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Abstract 

An  elastically  deformable  pressure  vessel  is  used  as  an  actuator. 

The  internal  pressure  is  controlled  hydraulically  via  force  applied  to  an 
external  bellows.  The  resulting  elongation  of  the  vessel  is  a  linear  function 
of  the  input  force,  and  depends  on  the  physical  properties  of  the  vessel, 
fluid,  and  bellows.  The  design  includes  an  orifice  through  which  fluid 
flows,  adding  damping  to  the  actuator. 

Mathematical  models  of  the  actuator  are  developed  which  relate 
performance  of  the  actuator  to  its  geometry  and  material  properties.  A 
static  analysis  yields  the  linear  relationship  between  the  commanded  input 
force  and  the  resulting  elongation  of  the  vessel.  A  model  of  the  passive  (no 
input  force)  response  of  the  actuator  indicates  it  will  act  as  a  passive 
damper.  The  active  response  (application  of  a  dynamic  input  force  yielding 
a  dynamic  elongation)  is  limited  in  frequency  by  the  orifice  damping. 
Strategies  for  optimizing  the  actuator  size  and  material  properties  are 
developed  from  the  models. 

Actuators  manufactured  with  differing  materials  and  fluids  were 
tested.  The  static  results  are  linear  and  matched  the  analysis.  The  passive 
results  and  the  model  predictions  confirm  that  the  configuration  of  the 
actuators,  as  built,  is  poor  for  passive  damping.  The  correlation  between 
the  model  and  active  experiments  is  excellent.  Actuation  bandwidth  is 
shown  to  be  selectable  by  selecting  orifice  size.  Tailoring  the  material 
properties  of  the  vessel  by  the  use  of  optimally  designed  composite 
laminates  results  in  a  factor  of  two  improvement  in  the  performance. 

The  actuator  is  useful  for  structural  control.  It's  performance  is 
comparable  to  other  available  actuators,  it  is  constructed  of  off-the-shelf 
hardware,  and  it  has  the  advantages  of  a  built-in  frequency  limit  and  easily 
customizable  performance  characteristics. 

Thesis  Advisor:  Professor  Hugh  L.  McManus 
Title:  Boeing  Assistant  Professor 


2 


Acknowledgement 

i 

In  the  short  two  years  that  I  have  been  here  (which  recently  have 
seemed  like  two  decades),  I  have  had  the  privelege  of  meeting  and  working 
with  so  many  brillant  good-natured  people. 

The  top  of  the  list  begins  with  my  advisor.  Professor  Hugh  McManus. 
Without  his  guidance,  his  time,  and  a  couple  of  cups  of  coffee  (that  he 
bought),  this  thesis  would  not  be  possible.  His  insight  and  ability  to  solve 
problems  were  a  charachteristic  that  I  admired  (and  hoped  a  little  has  rubbed 
off  on  me).  I  am  very  grateful.  Thank  you. 

To  everyone  at  the  Space  Engineering  Center  whom  I  have  gotten  to 
know  and  joked  with  thanks  for  keeping  me  sane.  The  help  of  so  many  of  you 
was  much  appreciated.  I  certanly  will  miss  the  lab.  (I  think,  thanks  to  the 
Air  Force,  I  will  reach  orbital  velocity  the  first  time  out.) 

To  Professor  Jack  Kerrebrock,  thank  you  for  the  concept  and  for 
recruiting  me.  I  hope  this  thesis  represents  just  the  first  tidbit  of  what  is 
possible  with  your  idea. 

To  my  wife,  who  married  me  knowing  the  thesis  was  on  the  horizon 
but,  probably,  did  not  anticipate  the  role  of  late-night  editor,  thank  you.  If  it 
had  not  been  for  you  my  organization  may  have  drowned  me.  I  consider 
myself  extremely  fortunate  to  have  found  a  person  as  understanding  and 
willing  to  help  as  you.  You  have  demonstrated  so  clearly  that  our  friendship 
and  companionship,  no  matter  what  the  future  entails,  is  a  constant  that  we 
can  trust  in.  Thank  you. 

This  research  was  funded  by  the  Space  Engineering  Research  Center 
(SERC)  at  the  Massachusetts  Institute  of  Technology.  I  am  grateful  to  have 
had  the  opportunity  to  work  on  this  project. 

I  also  appreciate  the  support  of  the  Air  Force  Institute  of  Technology, 
Wright  Patterson  Air  Force  Base,  Ohio. 


3 


Table  of  Contents 


Abstract .  2 

Acknowledgments .  3 

Nomenclature .  11 

Chapter  l:  A  Fluid  Elastic  actuator .  14 

1.1  Introduction .  14 

1.2  Actuator  Concept .  16 

1.3  Static,  Passive,  and  Active  Modes  of  Operations .  18 

1.4  Approach .  19 

1.5  Thesis  Outline .  19 

Chapter  2:  Structures  and  Actuators  in  Structures . 22 

2.1  Chapter  Outline .  22 

2.2  The  Controlled  Structures  Field .  22 

2.2.1  A  1-DOF  Model .  24 

2.2.2  A  Control  Approach .  27 

2.2.3  Passive  Damping  Approach .  28 

2.2.4  Active  Control  Approach .  30 

2.3  Previous  Applications  of  Pressure  Actuators  in  Structures .  34 

2.4  The  P-Strut  vs.  Previous  Pressure  Actuators .  35 

Chapter  3:  Developing  the  Performance  Equations . 36 

3.1  Chapter  Outline .  36 

3.2  Definition  of  the  P-Strut  Geometry .  38 

3.3  Derivation  of  the  Static  Equations .  40 

3.3.1  Constitutive  Equations .  40 

3.3.2  Pressure-Strain  Relationships .  46 

3.3.3  Derivation  of  the  Actuation  Authority,  Input  Force,  and 

Stroke  Requirements .  48 

3.3.4  Volumetric  Changes  and  Input  Force/Stroke 

Requirements .  51 

3.3.5  Static  Equations  for  an  Isotropic  Circular  Cylinder .  52 

3.3.6  Static  Equations  for  an  Anisotropic  Cylinder .  56 


4 


3.4  Analysis  of  the  Viscous  Fluid-Orifice  Damping .  58 

3.4.1  Laminar  Hagen-Poiseuille  Channel  Flow .  59 

3.4.2  Turbulent  Fully-Developed  Channel  Flow .  61 

3.4.3  Accelerating  and  Oscillating  Channel  Flow .  61 

3.4.4  Other  Fluid/Orifice  Considerations .  63 

3.5  Passive  Performance .  64 

3.5.1  An  Analysis  of  the  Mechanical  Behavior .  64 

3.5.2  Discrete  Passive  Stiffness  Model .  09 

3.5.3  Passive  Frequency  Characteristics .  71 

3.5.4  Parameteric  Study  of  the  Passive  Performance .  75 

3.5.5  Comparison  With  the  D-Strut .  77 

3.5.6  Comparison  With  Piezoelectric  Resistive-Shunted 

Struts .  78 

3.5.7  Revised  Passive  Model .  79 

3.6  Derivation  of  Active  Performance  Equations .  82 

3.6.1  Mechanical  Analysis  of  the  Active  Performance .  82 

3.6.2  Active  Frequency  Characteristics .  85 

3.6.3  Parameteric  Study  of  the  Active  Performance .  89 

3.7  Closing  Comments  on  the  Performance  Equations .  91 

Chapter  4:  Optimizing  the  P-Strut . 92 

4.1  Chapter  Outline .  92 

4.2  Optimizing  the  Physical  Dimensions .  92 

4.3  Optimizing  the  Material  Properties  with  Composites .  94 

4.4  Viscous  Fluid  and  Orifice .  97 

4.5  Hardware  Selection .  98 

4.6  Final  P-Strut  Design .  99 

4.7  “Sub-Optimal”  Aspects  of  the  P-Strut  Design . 103 

Chapter  5:  Manufacturing  and  experimental  Procedures . 105 

5.1  Chapter  Outline . 105 

5.2  Manufacturing  Process . 105 

5.2.1  Chemical  Milling  Process  for  Thin  Aluminum 

Cylinders . 105 

5.2.2  Composite  Cylinders  for  the  Hybrid  and  All-Composite 

P-Strut . 107 

5.2.3  Attaching  Endcaps . 110 

5.2.4  Filling  the  Cylinder  with  Viscous  Fluid . 110 


5 


5.3  Instrumentation  and  Experimental  Procedures . Ill 

5.3.1  Instrumentation . 112 

5.3.2  Static  Experimental  Procedure . 114 

5.3.3  Passive  Experimental  Procedure . 115 

5.3.4  Active  Experimental  Procedure . 116 

Chapter  6:  Experimental  Results  and  correlation . 117 

6.1  Chapter  Outline . 117 

6.2  Static  Performance  Evaluation . 118 

6.2.1  Static  Experimental  Results . 118 

6.2.2  Correlation  of  the  Static  Results . 119 

6.3  Passive  Performance  Evaluation . 129 

6.3.1  Passive  Experimental  Results . 129 

6.3.2  Correlation  of  the  Passive  Results . 129 

6.4  Active  Performance  Evaluation . 134 

6.4.1  Active  Experimental  Results . 134 

6.4.2  Correlation  of  the  Active  Results . 135 

6.5  Discussion  of  Results . 139 

Chapter  7:  Overview  and  Conclusions . 141 

7.1  Overview . 141 

7.2  Conclusions . 142 

7.3  Future  Work . 143 

References . 144 

Appendix  A:  End  Effects  and  Shell  bending  theory . 149 

A.l  The  Boundary  Effects  for  Isotropic  Circular  Cylinders . 149 

A.2  The  Boundary  Effects  for  Composite  Circular  Cylinders . 152 


6 


List  of  Figures 


1.1  P-Strut  Concept:  An  Elastic,  Axially  Deformable  Body 

Under  Internal  Fluid  Pressure .  17 

1.2  Static,  Active,  and  Passive  Modes  of  Operation .  17 

2.1  An  Overview  of  the  Controlled  Structures  Field .  23 

2.2  1-DOF  Model .  26 

2.3  Bode  Diagram  of  1-DOF  Model .  26 

2.4  A  Control  System  Diagram .  2B 

2.5  Passively  Damped  System  Model .  29 

2.6  System  Actuator  Model .  32 

3.1  Description  of  the  P-Strut  Geometry .  39 

3.2  P-Strut  Bellows  and  Orifice  Geometry .  41 

3.3  Definition  of  Longitudinal  (1)  and  Hoop  (2)  Coordinate 

Directions .  41 

3.4  Typical  Composite  Laminate  Layup .  44 

3.5  Pressurized  Circular  Cylinder .  47 

3.6  P-Strut  Radial  Swelling  and  End  Effects .  50 

3.7  Pressurized  Isotropic  Circular  Cylinder  End  Effects .  50 

3.8  Longitudinal  and  Radial  Volume  Changes .  53 

3.9  Total  Volume  Change  and  Input  Stroke .  53 

3.10  A  Circular  Orifice  with  Viscous  Fluid  Flow  and  the  Mechanical 

Dashpot  Representation .  60 

3.11  P-Strut  Schematic  in  Passive  Mode  of  Operation .  65 

3.12  Passive  Discrete  Stiffness  Model .  70 

3.13  Passive  Performance  ( Kpau )  Bode  Diagrams:  Magnitude  and 

Phase  Versus  Frequency .  73 

3.14  A  Schematic  of  the  D-Strut  .  7B 

3.15  Revised  Passive  Discrete  Stiffness  Model .  81 

3.16  Revised  Passive  Performance  Bode  Plots .  81 

3.17  P-Strut  Schematic  in  Active  Mode  of  Operation .  83 

3.18  Discrete  Active  Performance  (Admittance)  Model .  88 

3.19  Active  Performance  Admittance  (AJ  Bode  Diagrams: 

Magnitude  and  Phase  Versus  Frequency .  88 


7 


4.1  P-Strut  Vessel  Radius  Versus  Input  Force .  93 

4.2  P-Strut  Vessel  Length  Versus  Input  Force .  93 

4.3  Poisson’s  “Scissoring”  Effect:  Fc  Versus  Composite  Four 

Ply  Layups .  97 

4.4  P-Strut  Side  View . 100 

4.5  P-Strut  Axial  View . 101 

4.6  P-Strut  Top  View . 101 

4.7  Solenoid — Bellows — Orifice  Connections . 102 

4.8  An  “Optimal”  Vessel  for  the  P-Strut  Design . 104 

5.1  Chemical  Milling  Sample  Test  Results . 107 

5.2  Overlap  of  Composite  Plies  During  Layup  Cross-Sectional  View . 109 

5.3  Axial  Cross-Section:  P-Strut  Cure  Preparation . 109 

5.4  Bleed/Fill  Diagram . Ill 

5.5  P-Strut  Instrumentation . 113 

5.6  View  of  Component  Tester . 113 

6.1  Longitudinal  Strain  vs.  Hoop  Strain . 121 

6.2  Hoop  Strain  vs.  Input  Force . 121 

6.3  Pressure  vs.  Input  Force . 122 

6.4  Input  Stroke  vs.  Input  Force . 122 

6.5  Isotropic  Performance:  Elongation  vs.  Input  Force . 123 

6.6  Isotropic  Performance:  Elongation  vs.  Input  Stroke . 123 

6.7  Hybrid  Performance  with  10K  cs  Silicon  Fluid:  Elongation  vs. 

Input  Force . 124 

6.8  Hybrid  Performance  with  10K  cs  Silicon  Fluid:  Elongation  vs. 

Input  Stroke . 124 

6.9  Hybrid  Performance  with  30K  cs  Silicon  Fluid:  Elongation  vs. 

Input  Force . 125 

6.10  Hybrid  Performance  with  30K  cs  Silicon  Fluid:  Elongation  vs. 

Input  Stroke . 125 

6.11  Hybrid  Performance  with  Glycerol  Fluid:  Elongation  vs. 

Input  Force . 126 

6.12  Hybrid  Performance  with  Glycerol  Fluid:  Elongation  vs. 

Input  Stroke . 126 

6.13  Composite  Performance:  Elongation  vs.  Input  Force . 127 

6.14  Composite  Performance:  Elongation  vs.  Input  Stroke . 127 

6.15  Isotropic  Passive  Performance . 131 


8 


6.16  Hybrid  Passive  Performance . 131 

6.17  Hybrid  Passive  Performance  with  No  Fluid  . 132 

6.18  Composite  Passive  Performance . 132 

6.19  Isotropic  Active  Performance . 136 

6.20  Hybrid  Active  Performance  with  30K  cs  Silicon  and 

Glycerol  Fluids . 136 

6.21  Composite  Active  Performance  with  Increasing  Orifice 

Sizes  (30,  40,  55,  70,  82  mil) . 137 

6.22  Composite  Active  Performance  . 137 

A.l  Pressurized  Isotropic  Circular  Cylinder . 151 

A. 2  Hoop  and  Longitudinal  Strain  in  an  Isotropic  Cylinder . 153 

A. 3  Hoop  and  Longitudinal  Strain  in  a  Composite  Cylinder . 153 


9 


List  of  Tables 


3.1  Passive  Frequency  Characteristics .  72 

3.2  A  Parameteric  Study  of  the  Passive  Performance  Model .  76 

3.3  A  Parameteric  Study  of  the  Active  Performance  Model .  90 

4.1  A  Comparison  of  the  Numbers  of  Composite  Wraps  Versus 

The  Static  Performance  of  a  Hybrid  Aluminum-Composite .  96 

4.2  A  Comparison  of  P-Strut  Composite  Layups  Versus  the  Static 

Performance .  96 

4.3  Properties  of  the  P-Strut’s  Viscous  Fluids .  98 

4.4  The  Dimensions  of  the  Three  P-Strut  Designs . 103 

5.1  Instrumentation  List  and  Characteristics . 114 

6.1  Static  Performance  Results: . 128 

6.2  Constants  Used  in  Correlation  of  Static  Performance  Results . 128 

6.3  Measured  Passive  Performance  . 133 

6.4  Epoxy  Parameters  Fit  to  Passive  Data . 133 

6.5  Comer  Frequency  and  Active  Performance  Admittance . 138 


10 


Nomenclature 


A,X,Bd 
A 
A„ 


act 


'DC  act 


Bf 

B„ 


Cy,C„  K^.N.Y.e.r 


C, 

Clam 

Cfurb 

CL 

Do 

E 

E,  v,  G 

‘-'n  >  *-“22  >  v12  >  V21  >  “a 


1  L 

^iam’^turb 

h. 


I,tnd 


K 


State  Space  System  Matrices 

Composite  Thickness  Weighted  Stiffness  Matrix 

Admittance  of  P-Strut  in  Active  Model 

Admittance  at  DC  Frequency 

Area  of  Cylinder  Wall  (i.e.  Annulus  Area) 

Cross-Sectional  Area  of  Cylinder 

Effective  Cross-Sectional  Area  of  Bellows 

Cross-Sectional  Area  of  Orifice 

Bulk  Modulus  of  the  Fluid 

Bandwidth  (Frequency) 

Damping  Coefficient  Multiplied  by  Area  Ratio 
Control  Stete  Space  Matrices 
Damping  Coefficient 
Damping  of  Structure  or  Orifice  Modeled 
Damping  Coefficient 

Damping  Coefficient  Used  in  Passive  Model 

Laminar  Flow  Damping  or  Dashpot  Coefficient 

Turbulent  Flow  Damping  or  Dashpot  Coefficient 

Laminate  Compliance  Matrix 

Diameter  of  Orifice 

Material  Property  Matrix 

Isotropic  Material  Properties  (Young’s 

Modulus,  Poisson’s  Ratio,  Shear  Modulus) 
Laminate  Effective  Direction  Material  Properties 
Input  Force  into  Bellows 

Disturbance  Force  from  Structure  into  Actuator 
Force  from  Actuator  into  Structure 
Laminar  and  Turbulent  Head  Loss  Coefficients 
Head  Loss  Coefficient  Due  to  Roughness  or  Blunt 
Orifice  Ends 
Structural  Stiffness 


11 


k„k2 


K„ 


Kf 

Kl 


K 


POMM 


K 


past 


Ka 


Ks 

K, 


DCp 


K DCtp  >  Ke0  ,C, 


Kc 


ep’^ep 


cop 

L,  L 


fir 


Jdtu 


"o 

M 

P 


Ri,Re,Ra 


Re 


s 


t 

^cai  > 
^dtu 


t 


comp 


K 

Vo 

Z 


Stiffnesses  Associated  with  the  Passive  Damping 
Model 

Static  Stiffness  of  Bellows 

Static  Stiffness  of  Bellows  Multiplied  by  Area  Ratio 
Volumetric  Stiffness  of  Fluid 
Longitudinal  Stiffness 

Complex  Stiffness  of  P-Strut  in  Passive  Model 

Revised  Complex  Stiffness  of  Passive  Model 

Hoop  or  Radial  Dilation  Stiffness 

Structural  Stiffness  of  Bellows 

Passive  Stiffness  at  DC  Frequency 

Epoxy  Stiffness  and  Dashpot  Parameters 

Passive  Stiffness  at  Infinite  Frequency 

Physical  and  Effective  Length  of  Cylinder 

Entrance  or  Development  Length  in  Orifice 

Length  of  Orifice 

Mass  of  Structure 

Internal  Cylinder  Pressure 

Internal  Bellows  Pressure 

Interior  Radius  of  Cylinder,  Mean  Radius  of 

Cylinder,  and  Outer  Radius  of  Cylinder 

Reynold’s  Number 

Complex  or  Bode  Frequency  (ico) 

Thickness  of  Cylinder  Wall 

Thickness  of  Aluminum  and  Composite  Cylinders 
Time  Required  to  Develop  Laminar  Flow  in  Orifice 
Non-Dimensional  £*„ 

Time  to  Rise  (Step  Input) 

Velocity  of  Fluid  in  the  Orifice 
State  Space  Performance  Variable 


/?  Ratio  or  Frequency  of  Pole  Divided  By  Zero. 

Input  Bellows  Stroke 

SL  Elongation  or  Actuation  Authority 

Sr  Radial  Dilation  (Expansion) 


12 


ei  I»  &22 

Longitudinal  and  Hoop  Strain 

e 

Strain  Vector 

0 

Bode  Phase 

ri,  rf 

Loss  Factor  and  Peak  Loss  Factor 

<P 

Fiber  Orientation  Angle  for  Composites 

P 

Fluid  Absolute  Viscosity 

P 

Fluid  Density 

&IV  °22 

Longitudinal  and  Hoop  Stress 

a 

Stress  Vector 

u 

Fluid  Kinematic  Viscosity 

tn 

End  Effect 

CO 

Frequency 

v'd** 

Non-Dimensional  Oscillating  Fluid  Frequency  For 
Quasi-Laminar  Fluid  Flow 

of 

Peak  Damping  Frequency 

<4, 

Natural  Frequency 

Comer  Frequency 

<»p 

Zero  and  Pole  Frequencies 

c 

System  Damping  Coefficient 

AP 

Pressure  Loss  Through  Orifice  Shearing 

AVf 

Volume  Change  Due  to  Fluid  Compression 

Volume  Change  Due  to  Longitudinal  Expansion 

AVr 

Volume  Change  in  Radial  Direction  (Dilation) 

^ total 

Total  Volume  Change 

r 

Strain  Ratio:  Hoop  Strain  to  Longitudinal  Strain 

n,rc 

Isotropic  and  Composite  Strain  Ratios 

% 

Area  of  Ratio:  Cylinder  Area  to  Bellows  Area 

*0 

Area  of  Ratio:  Cylinder  Area  to  Orifice  Area 

Area  of  Ratio:  Bellows  Area  to  Orifice  Area 
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Chapter  1: 

A  Fluid  Elastic  Actuator 


1.1  INTRODUCTION 

The  established  trend  in  structural  design  is  towards  light  weight 
high  performance  structures.  Often  stringent  performance  criterion 
conflict  with  the  flexibility  inherent  in  light  weight  structures.  Flexible 
space  structures  under  current  development,  such  as  large  deployable 
reflectors,  optical  interferometers,  and  orbital  observation  platforms, 
exhibit  degraded  performance  in  the  presence  of  disturbance  vibrations. 

On  Earth,  the  engineering  fields  of  robotics,  high  performance  vehicles, 
and  precision  instruments  include  flexible  structures  which  require 
dynamic  stability  and  accurate  positioning  while  subjected  to  disturbance 
vibrations. 

An  adaptive  structure  employs  sensors  to  measure  performance  and 
actuators  to  correct  deviations.  Thus,  an  adaptive  structure  is  capable  of 
recognizing  and  manipulating  its  characteristics  or  states  [52].  In  order  to 
alter  the  states,  the  structure  must  possess  actuators  capable  of 
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commanding  inputs,  usually  displacement  or  force,  which  can  produce  the 
desired  and  corrective  action. 

Passive  devices  intrinsically  and  independently  react  to  their  local 
state.  A  damper,  for  example,  dissipates  energy  whenever  it  under  goes  a 
displacement.  Passive  devices  act  regardless  of  the  cause  of  the  local  state, 
and  without  knowledge  of  the  global,  structural  conditions.  The 
performance  of  a  passive  component  is  not  sensitive  to  desired  changes  in 
the  structure  and  cannot  be  optimized  or  controlled  once  installed. 

On  the  other  hand,  sensors  and  actuators  managed  by  an  automated 
control  system  can  accomplish  dramatic  improvements  in  performance. 
Errors  measured  by  sensors  and  feedback  to  the  control  system  may  be 
canceled  by  actuators.  In  addition,  the  physical  attributes  of  the  structure 
can  be  altered  (i.e.  a  part  could  be  stiffened  or  deployed/retracted)  via  an 
actuator  under  the  direction  of  a  control  system  [34]. 

Typically,  an  adaptive  structure  would  compliment  an 
actuator/sensor  control  system  with  passive  devices  in  order  to  achieve  the 
overall  performance  objective.  The  passive  component  reduces  the 
disturbance  and  the  level  of  control  effort  (i.e.  amount  of  actuated 
displacement  or  force)  needed  and  therefore  required  of  an  actuator. 

A  new  concept  of  commanding  displacements  or  force,  by  controlling 
the  deformation  of  a  fluid  containing  vessel,  results  in  an  actuator  with 
advantageous  characteristics.  This  idea  combines  an  active  component 
with  passive  fluid  damping. 
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1.2  ACTUATOR  CONCEPT 

A  closed  vessel  containing  a  pressurized  fluid  will  deform.  If  the 
vessel  is  elastic  and  the  fluid  pressure  is  controlled,  the  deformation  is 
elastic  and  predictable.  If  the  vessel  is  constrained  not  to  deform,  the 
pressure  response  is  a  force.  This  elastic  deformation  or  force  can  be  used 
as  an  actuator  to  command  displacement  or  force  into  a  structure.  The 
physical  and  mathematical  principles  which  explain  the  elastic 
deformation  of  a  pressurized  vessel  are  well  known.  Utilizing  a  fluid  as  the 
pressurized  medium  causes  the  dynamic  actuation  to  be  a  function  of  the 
viscous  fluid  motion,  as  well  as  the  vessel  properties.  This  viscous  fluid 
motion  produces  innate  passive  damping. 

Figure  1.1  depicts  the  geometry  of  a  pressurized  vessel.  The  internal 
pressure  causes  an  elastic  deformation,  SL,  or  actuation  force,  FL,  which 
depend  on  the  properties  and  dimensions  of  the  vessel.  An  input  force,  Fb, 
controls  the  pressure  via  a  piston-like  mechanism.  The  fluid  subjected  to  Fb 
is  connected  to  the  fluid  in  the  vessel  through  a  restricted  passageway  or 
orifice. 

The  actuation  displacement,  SL,  is  a  linear  function  of  Fb.  If 
constrained,  an  actuation  force,  FL,  is  a  linear  function  of  Fb  multiplied  by  a 
ratio  of  the  areas  and  dependent  on  the  vessel  properties.  This 
amplification  effect  is  analogous  to  a  hydraulic  jack  or  a  mechanical  lever 
where  a  large  pressure  force  is  gained  at  the  expense  of  a  large  input 
displacement.  As  the  viscous  fluid  squeezes  through  the  orifice,  the  fluid 
undergoes  a  shearing  action  and  releases  energy.  This  loss  of  energy  is  the 
source  of  damping. 
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FIGURE  1.1:  P-Strut  Concept:  An  Elastic,  Axially  Deformable 
Body  Under  Internal  Fluid  Pressure 


FIGURE  1.2:  Static,  Active,  and  Passive  Modes 
Of  Operation 
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By  tailoring  the  material  properties  of  the  vessel,  the  deformation  or 
the  actuation  force  can  be  optimized.  Furthermore,  the  type  of  actuation — 
elongation,  bending,  or  torsion — are  functions  of  the  geometry  and  material 
properties  of  the  actuator. 


1.3  STATIC,  PASSIVE,  AND  ACTIVE  MODES  OF  OPERATIONS. 

The  pressure  actuator  operates  in  three  distinct  modes  which  are 
illustrated  in  Figure  1.2.  During  static  operation,  where  the  velocity  of  the 
viscous  fluid  is  negligible,  the  elongation,  SL,  is  linear  with  the  input  force, 
Fb.  Thus,  the  actuator  could  command  and  hold  precise  structural  inputs 
by  maintaining  a  constant  fluid  pressure. 

In  the  presence  of  dynamic  structural  disturbances,  Fd,  the  actuator 
would  be  expanded  or  compressed  generating  SL  independent  of  Fb.  In 
order  to  accommodate  the  volume  changes  resulting  from  the  disturbance, 
the  fluid  in  the  vessel  would  be  pumped  or  sucked  through  the  orifice.  The 
fluid  shearing  in  the  orifice  would  cause  a  damping  of  the  disturbance. 

This  mode  of  operation  is  passive. 

The  third  mode  of  operation  is  active.  If  the  input  force,  Fb,  is 
dynamic,  the  fluid  in  the  reservoir  would  be  forced  into  the  orifice  at  a  rate 
directly  dependent  on  Fb .  However,  the  viscous  fluid  is  inhibited  by  the 
restricted  passageway.  Therefore,  Fb  applied  slowly  would  allow  the  fluid  to 
pass  through  the  orifice  and  into  the  vessel  (i.e.  near  static  conditions);  but 
with  a  rapidly  applied  Fb  the  fluid  would  not  squeeze  through.  Thus,  the 
actuator  cannot  command  displacement  or  force  at  high  rates  or 
frequencies.  This  physical  limitation  prevents  the  actuator  from  exciting 
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the  higher  dynamics  of  a  structure  or  amplifying  back  into  a  structure 
control  system  noise. 

The  technology  needed  to  create  and  implement  the  pressurized  fluid 
elastic  actuator  is  available  from  “off-the-shelf"  hardware.  The  concept  is 
straightforward  and  an  inexpensive,  simple  alternative  to  other  actuator 
approaches. 


1.4  APPROACH 

The  proposed  pressurized  fluid  elastic  actuator  (dubbed  the  P-Strut) 
was  designed  with  the  objective  of  achieving  actuation  comparable  to  other 
available  actuators.  The  dimensions  of  the  P-Strut  were  optimized  by 
solving  the  pressure-displacement  relationships  for  pressure  vessels.  The 
investigation  examined  the  advantages  of  tailoring  the  material  properties 
to  maximize  the  axial  elongation.  Analytical  models  were  developed  to 
represent  the  three  modes  of  operation. 

Three  P-Struts  were  manufactured  to  prove  the  concept  and 
demonstrate  the  benefits  of  customizing  the  material  properties.  The 
actuators  were  tested  in  each  of  the  three  modes  operation  for  a  variety  of 
fluids  and  orifice  diameters.  The  results  correlated  well  with  the  models. 


1.5  THESIS  OUTLINE 

This  thesis  describes  the  development  and  verification  of  a  fluid 
elastic  actuator. 
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Chapter  2  provides  the  background  to  controlled  structures  and 
discusses  pneumatic  piston  devices  which  have  been  previously  applied  to 
the  control  of  flexible  structures. 

The  analytical  groundwork  for  the  design  is  presented  in  Chapter  3. 
The  investigation  begins  by  developing  the  pressure-displacement 
relationships.  The  static  performance  equations,  defined  as  input  force,  Fb, 
and  input  stroke,  Sb,  versus  actuator  displacement,  SL,  are  derived  in  terms 
of  the  P-Strut's  material  properties  and  physical  dimensions.  The  viscous 
fluid  and  the  orifice  interaction  is  considered  from  solutions  to  channel  flow 
problems.  The  three  modes  of  operation  are  described  via  models  which 
mimic  the  physics  of  the  P-Strut.  The  parameters  of  these  models  include 
elastic  material  and  volumetric  fluid  stiffnesses,  as  well  as  a  rate 
dependent  parameter  (i.e.  a  dashpot)  for  the  fluid-orifice  effect.  The  models 
are  used  to  define  the  appropriate  mathematical  equations.  A  parameteric 
study  of  the  models  emphasize  the  important  parameters.  The  passive 
model  was  updated  to  account  for  an  additional  damping  source  discovered 
during  testing. 

The  equations  developed  in  Chapter  3  are  combined  to  determine  the 
optimal  actuator  dimensions  in  Chapter  4.  An  isotropic,  a  composite- 
aluminum  hybrid,  and  an  all-composite  P-Strut  are  described  and 
optimally  designed  to  demonstrate  the  advantages  of  tailoring  the 
actuator’s  material  properties.  The  “ofF-the-shelf”  components  of  the  P- 
Strut  are  discussed  as  well  as  the  properties  of  the  viscous  fluids  utilized. 

Chapter  5  details  the  manufacturing  procedures  and  the  test 
methods.  The  construction  included  the  chemical  milling  of  aluminum 
tubes,  and  the  layup  of  aluminum-composite  (hybrid)  and  all-composite 
tubes.  The  method  of  filling  the  actuator  with  fluid  is  presented.  The  test 
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objectives  and  setup  are  described.  The  techniques  and  equipment  utilized 
to  determine  the  static,  passive,  and  active  performance  of  the  P-Strut  are 
discussed. 

The  experimental  results  are  reported  and  correlated  in  Chapter  6. 
The  performance  of  the  three  P-Struts  with  a  variety  of  orifice  and  fluid 
combinations  is  presented.  The  static  performance  was  linear.  The  active 
results  compared  favorably  with  the  static  results  and  the  model  developed 
in  Chapter  3.  The  advantage  of  tailoring  the  material  properties  is 
supported  by  the  data.  Passive  damping,  probably  caused  by  the  epoxy  used 
in  manufacturing,  was  discovered  during  testing.  This  damping 
overshadowed  the  orifice  damping. 

A  summary  of  the  scientific  and  engineering  benefits  of  this  study 
conclude  the  thesis  in  Chapter  7. 
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Chapter  2: 

Structures  and  Actuators  in  Structures 


2.1  CHAPTER  OUTLINE 

This  chapter  presents  a  short  review  of  flexible  structures  and  the 
relevant  aspects  of  structural  control.  A  one  degree-of-freedom 
(1-DOF)  model  is  used  to  explain  issues  involved  in  passive  damping  and 
active  control.  The  use  of  an  actuator  in  a  feedback  control  system  is 
examined.  Three  general  categories  of  feedback  control  are  considered. 

Previous  actuators  which  have  used  pressure  as  an  actuation  source 
are  discussed  and  compared  with  the  P-Strut  concept.  The  differences 
between  a  piston  pneumatic  actuator  and  the  fluid  elastic  actuator  are 
clarified.  Structural  control  issues  surrounding  the  proposed  fluid  elastic 
actuator  are  presented. 

2J2  THE  CONTROLLED  STRUCTURES  FIELD 

Figure  2.1  depicts  an  overview  of  the  controlled  structures  field. 
Generally,  a  structure  is  first  modeled  with  a  set  of  assumptions  and 
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FIGURE  2.1:  An  Overview  of  the  Controlled  Structures  Field 
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estimated  physical  parameters.  In  reference  10,  Craig  published  a  broad 
survey  of  recent  modeling  topics  and  techniques.  Modeling  provides  means 
to  predict  the  performance  of  the  structure  and  design  a  control  approach. 

2.2.1  A 1-DOF I J3DEL 

A  1-DOF  system  consisting  of  a  single  mass,  spring,  and  dashpot  is 
pictured  in  Figure  2.2.  The  model  represents  a  structure  with  mass,  AT, 
elastic  stiffness,  K,  and  material  damping,  C,  subjected  to  a  disturbance 
force,  Fd.  This  system  is  useful  for  illustrating  several  significant 
structural  control  concepts.  The  mathematical  representation  of  the 
system  is  the  equation  of  motion: 

Mx  +  Cx  +  Kx  =  Fd  (2.1) 

or: 

F 

x  +  2Z(Onx  +  Q>l=-*r  (2.2) 


IK  C 

where:  <a„  =  ,  —  and  {  = - the  structure’s  natural  frequency 

V  M  2  (onM 

and  damping  coefficient,  respectively.  £  is  a  function  of  the  material 
damping  parameter,  C,  which  in  turn  is  dependent  on  the  material  type, 
and  in  some  instances,  on  the  loading  [11,41,48].  For  a  typical  structure,  £ 
is  very  small  (0.1-1.0%). 

An  analysis  of  the  model  is  usually  accomplished  in  the  frequency 


domain  or  the  time  domain  [37,51].  In  the  time  domain,  assuming  the 
damping  is  negligible  (£  =0),  the  equation  of  motion  can  be  written  as  two 
first  order  differential  equations  which  together  are  called  the  state  space 
equation: 
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X 


(2.3) 


j®,}^  =  (X)  =  [A](X)  +  [Ba]FJ 


X  is  the  state  vector  (i.e.  displacement  and  velocity),  A  is  the  system  matrix, 
and  Bd  is  the  disturbance  force  mapping  matrix  which  indicates  how  the 
force  is  directed  into  the  states. 

A  frequency  domain  representation  can  be  derived  by  taking  the 
Laplace  Transform  of  2.2  with  zero  initial  conditions: 


X(s) 

~Fd(s) 

K 


O). 


s2+0)2n 


=  G(s) 


(2.4) 


where  G(s )  relates  the  dynamic  displacement  X(s)  to  dynamic  input  Fd(s) 
and  the  system  parameters  [37,51]. 

G(s)  is  the  state  transfer  function  [37,51].  The  roots  of  the 
denominator  are  called  poles  and  the  roots  of  the  numerator  are  called 
zeroes.  Two  repeated  poles  are  at  (on .  The  physical  interpretation  of  a  pole 
and  zero  can  be  observed  from  Figure  2.3,  the  Bode  Diagram  of  G(s)  [37,51]. 
The  Bode  magnitude  and  phase  (i.e.  the  magnitude  and  phase  angle  of  the 
sinusoidal  response,  X,  to  a  sinusoidal  input  force,  Fd)  are  depicted  versus 
the  complex  frequency  s.  The  peak  in  the  magnitude  and  the  shift  in  the 
phase  occurs  at  con .  This  frequency  location  is  labeled  the  natural 
frequency,  resonant  frequency,  or  frequency  of  the  structural  mode.  This 
frequency  corresponds  to  a  particular  motion  or  mode  shape  of  the 
structure.  For  the  system  in  Figure  2.2  the  mode  shape  is  the  mass  moving 
left  to  right.  At  an  undamped  pole  the  response  to  a  disturbance  would  be 
infinite.  If  material  damping  is  included  (£?0),  the  height  of  the  peak  and 
the  sharpness  of  the  phase  shift  are  dependent  on  the  inverse  of  £  [12,37,51], 
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FIGURE  2.3:  Bode  Diagram  of  1-DOF  Model 


At  frequencies  approaching  infinity  the  response  rolls  off  to  zero  at  an 
average  slope  on  a  log-log  scale  of  approximately  the  number  of  poles  minus 
the  number  of  zeroes  [12]. 

2JL2  a  Control  approach 

The  purpose  of  a  control  system  is  usually  to  minimize  the 
performance  deviations  and  therefore  minimize  the  associated  deviations 
in  the  structural  states.  This  amounts  to  controlling  the  modes  of  the 
system  or  reducing  the  resonant  peaks  of  the  transfer  function.  Usually, 
the  approach  to  control  design  is  to  damp  or  move  the  resonant  peaks  which 
are  critical  to  the  performance. 

To  accomplish  this  objective  the  plant  or  modeled  system  must  be 
paired  with  a  control  approach.  The  scheme  could  be  as  simple  as  adding 
passive  components  or  designing  a  feedback  control  system  using  actuators 
and  sensors.  The  type  and  level  of  control  depends  on  the  needed 
performance. 

Figure  2.4  depicts  a  typical  feedback  control  diagram.  The  states  are 
labeled  X  and  are  observed  as  the  output  of  G(s )  caused  by  the  input  of  the 
combined  disturbance  force,  Fd,  and  an  actuator,  Fa.  Y  are  the  sensor 
outputs  which  are  related  to  the  states  of  the  structure  by  a  gain  matrix  Cy 
and  are  corrupted  by  a  sensor  noise  N.  The  performance  of  the  system,  Z,  is 
given  by  C.Xand  consists  of  those  states  which  should  be  controlled.  The 
compensator  loop  gain  and  the  dynamics  of  the  actuator  and  sensors  are 
combined  in  the  compensator/actuator  block,  K^is).  The  error,  e,  equals 
the  difference  between  Y  and  the  reference,  r.  If  the  gain  and  state 
matrices  are  linear  and  constant,  the  system  is  linear  and  time  invariant 
(abbreviated  LTI)  [37,51]. 
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FIGURE  2.4:  A  Control  System  Diagram 


To  reduce  the  performance  deviations  as  observed  by  Z,  several  types 
of  control  strategies  are  available,  and  research  continues  into  new  control 
strategies.  However,  for  purposes  of  this  discussion  four  generally  well 
known  approaches  are  examined:  passive  damping,  rigid  body  control,  low 
authority  control,  and  high  authority  control. 

2.2.3  passive  Damping  Approach 

Passive  damping  reduces  the  peak  of  the  performance  transfer 
function  by  increasing  the  system  damping,  £.  In  most  instances,  passive 
damping  is  added  before  feedback  control  and  is  thus  independent  of  the 
control  system.  Passive  damping  changes  the  system  transfer  function 
G(s). 

Figure  2.5  depicts  a  new  structural  model  where  a  generic  passive 
damping  device  has  replaced  the  structural  stiffness  or  spring  element  K  in 
Figure  2.2.  The  damper  is  modeled  as  a  pair  of  springs  and  a  dashpot. 
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This  model  of  a  damper  is  representative  of  the  passive  mode  of  operation 
for  the  P-Strut.  The  device  is  in  the  load  path  of  the  structure  and  therefore 
must  support  the  mass. 


Passive  Damper 

Figure  2.5:  Passively  Damped  System  Model 

The  purpose  of  the  damper  is  to  add  damping  to  the  structure  over  a 
particular  bandwidth.  The  dashpot  and  springs  have  a  characteristic 
damping  curve  which  offers  significant  damping  within  a  limited 
bandwidth.  Inserting  the  device  into  the  structure  increases  the  damping 
coefficient  of  the  system  within  the  bandwidth  of  the  device.  In  order  to 
achieve  the  optimal  reduction  in  the  resonant  peak,  the  frequency  of 
maximum  damping  within  the  bandwidth  should  be  designed  to  match  the 
resonant  frequency.  However,  if  the  resonant  frequency  moves  from  the 
peak  damping  frequency,  the  damper  is  mistuned.  The  result  is  less, 
possibly  insufficient,  damping. 

Damping  is  achieved  by  dissipating  the  energy  in  a  system.  Energy 
can  be  removed  from  the  system  through  a  strain  mechanism.  For 
example  a  constrained  viscoelastic  layer  attached  to  a  structural  member 
dissipates  energy  by  shear  strain.  This  mechanism  offers  reasonable 
damping  with  little  mass  penalty  [3,39,49].  The  Honeywell  D-Strut — a 
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dashpot  damper  currently  in  operation  aboard  the  Hubble  Space 
Telescope — provides  significant  damping  over  a  reasonable  band  of 
frequencies  by  displacing  a  fluid  through  a  restricted  passageway 
[2,9,14,35].  This  damping  concept  is  analogous  to  the  P-Strut.  A  Piezo* 
resistive-shunted  device,  which  when  strained  generates  a  voltage  that  is 
dissipated  in  a  resistor,  also  offers  damping  over  a  moderate  frequency 
range  [20].  These  types  of  passive  devices  have  a  fixed  frequency  at  which 
the  damping  is  maximum  and  a  relatively  broad  band  about  this  peak  in 
which  the  damping  is  effective.  In  contrast,  resonant  dampers  such  as 
proof  mass  dampers  or  piezo-inductive-resistive-shunted  devices  are 
capable  extracting  a  significant  amount  of  energy  and  can  be  tuned  to  a 
particular  frequency.  However,  these  devices  add  additional  modes  to  the 
system,  and  are  only  effective  in  a  narrow  band  about  the  tuned  frequency 
[20,29]. 

Although  passive  components  are  benign  (i.e.  they  cannot  add 
energy),  they  are  only  capable  of  reacting  to  locally  perceived  displacements 
or  forces  and  are  independent  of  and  unaware  of  the  global  performance. 

2.2.4  Active  Control  approach 

An  adaptive  structure  requires  control  systems  which  can  sense  and 
command  a  structural  change  and  compensate  for  unwanted  changes. 
Figure  2.6  depicts  the  1-DOF  Model  with  an  actuator  in  parallel  with  the 
structural  stiffness.  The  active  controller  or  compensator,  K^is),  in  Figure 
2.4  which  controls  the  feedback  gains  and  thereby  the  actuator  response,  F„ 
depends  on  the  type  or  level  of  control  needed. 
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2.2.4. 1  Rigid  Body  Control 

Rigid  body  control  involves  deploying,  retracting,  or  positioning  an 
object.  In  order  to  accomplish  a  rigid  body  movement,  an  input 
displacement  is  needed.  If  the  mass  is  independent  from  the  wall  in  the  1- 
DOF  model  (i.e.  /l=0  in  Figure  2.6),  a  commanded  actuator  displacement 
would  cause  the  mass  to  move  to  a  new  position  and  hold.  If  the  mass  were 
truly  rigid  and  no  disturbance  force  was  present  the  movement  would  be 
precise.  However,  if  K*Q,  or  the  mass  has  additional  flexibilites  not  shown 
in  the  figure,  the  movement  would  involve  a  decaying  vibration  or  ring 
down  around  the  displaced  position.  Preferably,  an  actuator  would  not  add 
energy  to  the  flexibility  of  the  system.  This  can  be  accomplished  either  by 
input  shaping,  which  diminishes  the  energy  transmitted  to  the  flexible 
modes  of  the  structure  by  negatively  reinforcing  the  resonant  frequencies 
[45],  or  by  an  actuator  which  rolls  off  at  frequencies  lower  than  those  of  the 
flexibility  of  the  mass.  In  a  feedback  control  system,  the  measured  states 
can  be  used  to  adjust  the  inputs  and  reposition  or  realign  the  structure  via 
the  actuator  especially  in  the  presence  of  disturbance  force,  Fd. 

Ideal  actuators  are  linear  such  that  a  feedback  measurement  yields  a 
linear  correction  through  the  compensator.  However,  piezoceramics 
exhibit  mild  hystersis  effects  and  electrostrictives  are  quadratic  (or  higher) 
with  respect  to  the  driving  voltage  [4,13,32,50]. 
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FIGURE  2.6:  System  Model  Actuator 


2J2.4.2  Low  Authority  Control 

Low  authority  control  (LAC)  decreases  the  resonant  peaks  in  Figure 
2.2  by  actively  adding  damping  to  the  structure.  The  objective  is  similar  to 
that  of  passive  damping;  however,  the  achievable  peak  damping  is  higher 
and  the  frequency  can  be  tuned.  A  proportional  derivative  (PD) 
compensator  feeds  back  a  measurement  of  the  displacement  and  the 
velocity  to  influence  both  the  natural  frequency  and  system  damping  [37,51], 
By  properly  designing  the  compensator,  the  structure  can  be  significantly 
damped.  A  linear  quadratic  regulator  (LQR),  which  is  a  full  state  feedback 
controller,  accomplishes  the  peak  damping  for  a  single  mode  system  in  the 
same  manner  as  the  PD  controller  [12].  LQR  is  one  type  of  modern  control 
design  technique.  An  advantage  of  modern  controllers  is  the  use  of  the 
state  space  representation  (equation  2.3)  which  is  optimal  for  computer 
implementation.  The  deviations  in  performance  can  be  written  as  a  cost 
which  can  be  minimized  by  mathematical  operations  resulting  in  a  stable 
compensator  design.  [12,37,51]. 

Regardless  of  the  compensator,  the  actuator  must  be  capable  of 
commanding  an  input  at  the  frequency  of  the  mode  which  is  to  be  damped. 
The  frequency  bandwidth  over  which  the  component  can  actuate  dictates  to 
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what  frequency  and  to  what  extent  active  damping  is  possible.  If  the 
actuator  cannot  command  force  or  displacement  at  a  particular  frequency 
no  damping  or  control  is  possible.  Therefore,  care  must  be  exercised  in 
selecting  the  appropriate  actuator.  Piezoceramic  materials  and 
electrostrictives  are  capable  of  excitation  over  a  broad  range  of  frequencies 
[46].  At  frequencies  out  to  IK  Hz,  Scribner,  et.al.  employed  a  piezoelectric 
component  to  isolate  transmitted  disturbances  in  a  model  of  a  rotor  blade 
[44].  At  lower  frequencies  (<  25  Hz),  Hallauer,  et.al.  used  an  air  jet  thruster 
alone  and  paired  with  a  reaction  mass  actuator  to  damp  a  planar  truss 
[22,23]. 

2J2.4.3  High  authority  Control 

High  authority  control  (HAC)  involves  shifting  the  structural 
frequencies  and  altering  the  mode  shapes  of  the  structure  [12].  HAC 
means  moving  the  peak  and  lowering  the  entire  curve  in  Figure  2.3 
therefore,  greatly  diminishing  the  performance  deviations.  The  magnitude 
of  the  commanded  input  must  be  comparable  to  the  disturbance.  In 
essence,  HAC  attempts  to  cancel  the  disturbance  such  that  from  Figure  2.6, 
Fa--Fd.  HAC  can  involve  a  multi-input-multi-output  problem  with  many 
actuators  and  sensors.  HAC  compensators  can  be  determined  from 
frequency  weighted  cost  functions  including  techniques  to  handle 
parametric  uncertainties  and  unmodeled  dynamics  [12]. 
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2 J3  PREVIOUS  APPLICATIONS  OF  PRESSURE  ACTUATORS  IN 


STRUCTURES 

Due  to  the  inherent  bandwidth  limitations  on  pneumatic  or  pressure- 
type  actuators,  structural  applications  have  been  limited.  Sievers  and 
VonFlotow  appropriately  dismissed  pneumatic  actuators  in  developing  an 
isolator  for  damping  acoustic  vibrations  [46].  Nevertheless,  pneumatic 
actuators  have  been  successfully  used  at  low  frequencies. 

Rafati  designed  and  tested  a  pneuma  piston  actuator  with 
application  to  passenger  trains  [40].  His  component  was  utilized  in  a  LAC 
system  to  damp  modeled  disturbances  transmitted  f  Dm  the  rail  to  the 
train.  These  disturbances  included  the  accelerations  due  to  train  rocking 
and  bouncing  and  occurred  at  frequencies  below  5  Hz.  The  opening  and 
closing  delays  in  a  solenoid  valve  limited  the  performance  of  his  system. 

The  pneumatic  air-jet  developed  by  Hallauer  and  Smith — used  to 
damp  a  planar  truss — expelled  pressurized  air  to  generate  a  force  (i.e. 
thrust).  The  air-jet  had  a  limited  bandwidth  due  to  delays  in  opening  and 
dosing  the  air  valve  [22]. 

Lim,  et.al.  used  unspecified  pneumatic  actuators  in  comparing  three 
multi-input-multi-output  control  systems.  The  structure  consisted  of  a 
fifty-one  foot  truss  connected  to  a  sixteen  foot  diameter  reflector  with  the 
highest  controlled  frequency  of  1.87  Hz  [30].  His  results  suggested  that  a 
PD  type  controller  with  feedback  to  the  pneumatic  actuators  compared 
favorably  to  other  controllers. 

Kamopp  proposed  an  active  and  passive-active  isolation  system  for  a 
high-speed  ground  vehicle  using  a  modem  control  theory  approach  [26]. 
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Butsuen  proposed  a  semi-active  isolation  system  to  isolate  roadway 
disturbances  without  severely  limiting  an  automobile's  handling 
characteristics  [7].  Since  a  substantial  force  is  required  over  a  modest 
bandwidth,  these  control  systems  could  involve  a  fluid  elastic  actuator. 


2.4  THE  P-STRUT  VS.  PREVIOUS  PRESSURE  ACTUATORS 

The  proposed  P-Strut  commands  a  force  or  displacement  by 
controlling  the  internal  pressure  in  an  elastic  deformable  cylinder.  On  the 
other  hand,  the  piston  actuator  developed  by  Rafati  used  a  solenoid  to  open 
and  close  a  pressurized  chamber.  Hallauer’s  air-jet  released  a  directed 
stream  of  high  pressure  air  to  generate  on/off  thrust. 

The  proposed  fluid  elastic  actuator  has  an  inherent  frequency  roll-off 
similar  to  other  pneumatic  actuators.  However,  in  the  P-Strut  the 
bandwidth  is  limited  by  the  fluid-orifice  interaction,  whereas  in  piston 
actuators  and  the  air-jet,  the  frequency  bandwidth  is  limited  by  the  opening 
and  closing  of  the  solenoid  valves  [40,22].  Moreover,  the  viscous  fluid  adds 
passive  damping  qualities.  When  not  actuating,  the  P-Strut  is  capable  of 
absorbing  energy,  thus  further  improving  the  performance  of  the  structure. 
In  contrast,  the  air-jet  and  piston  actuators  are  not  designed  to  add  passive 
damping. 
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Chapter  3: 

Developing  the  Performance  Equations 


3.1  CHAPTER  OUTLINE 

This  chapter  examines  the  physical  principles  of  the  P-Strut  and 
develops  the  mathematical  equations  which  describe  the  actuator's  three 
modes  of  operation:  static,  passive,  and  active. 

Section  3.2  describes  the  dimensions  of  the  P-Strut  that  are  relevant  to 
the  analyses. 

In  Section  3.3,  the  pressure-strain  relationships  are  developed  via  the 
stress-strain  constitutive  equations  for  isotropic  and  anisotropic  materials. 
The  axial  elongation  of  the  structure  is  determined  as  a  function  of  the 
input  force  and  input  stroke.  These  equations  define  the  static  performance 
of  the  actuator.  Anisotropic  material  properties  are  shown  to  enhance  the 
performance. 

Section  3.4  investigates  the  damping  interaction  between  the  viscous 
fluid  and  the  orifice.  An  effective  loss  coefficient  is  estimated  by 
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considering  laminar  and  turbulent  channel  flows.  In  addition, 
accelerating  and  sinusoidally  varying  pressure  flows  are  examined. 

A  mechanical  analysis  details  the  passive  mode  of  operation  and 
leads  to  a  three  parameter  model  in  Section  3.5.  The  parameters  include 
discrete  elastic  and  volumetric  stiffnesses,  as  well  as  a  dashpot 
representation  of  the  fluid-orifice  interaction.  The  physical  and 
mathematical  development  indicates  that  these  stiffnesses  compete  against 
one  another  reducing  the  passive  damping.  The  P-Strut  performance  is 
compared  to  the  Honeywell  D-Strut  and  piezoelectric  resistive-shunted 
struts.  The  passive  model  is  updated  to  include  the  additional  damping 
encountered  during  testing.  The  revised  model  based  on  the  measured 
performance  has  damping  which  overwhelms  the  predicted  orifice 
damping. 

Section  3.6  examines  the  active  performance  of  the  P-Strut.  A  second 
discrete  stiffness  model  is  developed  with  similar  physical  parameters. 
Mathematically  and  physically,  the  model  defines  the  non-collocated 
transfer  function  between  the  input  force  and  the  actuation  displacement. 
The  static  performance,  as  well  as  the  frequency  roll-off  due  to  the  fluid  and 
orifice,  are  properly  incorporated  into  the  model.  A  parameteric  study 
demonstrates  the  important  parameters. 
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3.2.  DEFINITION  OF  THE  P-STRUT  GEOMETRY 


Figure  3.1  depicts  a  schematic  of  the  P-Strut.  The  fluid-elastic 
actuator  is  a  circular  cylindrical  pressure  vessel  which  is  sealed  with  rigid 
flat-plate  endcaps.  The  fluid  inside  the  cylinder  is  connected  via  an  orifice 
to  the  fluid  inside  the  reservoir  or  flexible  bellows. 

The  cylindrical  actuator  has  a  length  L,  with  inner  radius  Rit  outer 
radius  R0,  and  thickness  t.  Re  represents  the  average  radius  or  (V2)(Ri  +Ra) 
Since  t«Re,  the  three  radii  are  approximately  equal  (i.e.  R0~Rc=Rl).  Thus, 
the  cross-sectional  area,  Ac,  and  the  annulus  or  vessel  wall  area,  Ax,  can 
be  expressed  as: 


Ae  =  *  nR2e 


Aac  =  x(R20-R*)~2jtR<t 


(3.1) 


The  bellows  is  depicted  in  Figure  3.2.  The  reservoir  or  bellows  has  an 
effective  cross-sectional  area  Ab.  This  dimension  is  difficult  to  compute 
from  the  geometry  of  the  bellows.  (In  this  study,  Ab  is  determined  directly 
from  the  experimental  results.) 

The  orifice  has  cross-sectional  area  A0  and  length  L0. 

The  ratio  of  the  P-Strut  cross-sectional  area,  Ac,  to  the  effective 
bellows  area,  Ab,  equals  the  fluid  lever  introduced  in  Chapter  1.  This  fluid 
lever  is  represented  by  *P6.  The  ratio  of  Ac  to  the  area  of  the  orifice,  A0,  is 
defined  as  *F„.  Ab  divided  by  A0  equals  *f'o6.  As  equations: 

•p.  =A  =  —  '?06=  —  =  —  (3.2) 

4  A,  A  'f'b  A 
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FIGURE  3.1:  Description  of  the  P-Strut  Geometry 
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Figure  3.3  depicts  the  element  ABCD  extracted  from  the  surface  of 
the  P-Strut  cylinder  illustrated  in  Figure  3.1.  The  longitudinal  direction  is 
labeled  L  or  1  and  the  hoop  direction  is  labeled  H  or  2.  The  curvature  of  the 
element  in  the  hoop  direction  is  Rc.  The  element  is  not  curved  with  respect 
to  the  longitudinal  direction. 

3^  DERIVATION  OF  THE  STATIC  EQUATIONS 

As  established  in  Chapter  1,  the  P-Strut  is  an  elastically  deformable 
pressure  vessel.  The  static  analysis  derives  the  elongation  of  the  P-Strut 
(i.e.  actuation  authority)  as  functions  of  the  input  force  and  stroke. 

3.3.1  Constitutive  Equations 

In  general,  the  application  of  an  arbitrary  load  induces  normal  and 
shear  stresses  in  the  cylinder  wall  which  can  be  described  in  the  1,  2,  and  3 
coordinate  directions.  However,  since  the  thickness  of  the  cylinder,  t ,  is 
much  smaller  than  the  radius,  Rc,  the  stresses  through  the  thickness  are 
negligible  (i.e.  al2,  o23 ,  and  <rM=0).  With  this  plane-stress  state  assumption, 
the  constitutive  equations  are: 

^1111  ^1122  2Eim  f£n 

'  ^22  ’  =  ^1122  ^2222  ^^22X2  j  ^22  '  (3.3) 

.^12,  .^1112  ^2212  2Ej2j2_  [^12  , 

or  abbreviated: 

(o}*[E]{e)  (3.4) 
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FIGURE  3.2:  P-Strut  Bellows  and  Orifice  Geometry 
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The  components  of  the  E  matrix  represent  the  cylinder’s  material 
properties.  The  2  appears  in  the  last  column  of  E  since  the  tensor  strain, 
e12,  instead  of  the  engineering  strain,  Y12,  is  used. 

If  a  pair  of  orthogonal  planes  of  material  property  symmetry  exists, 
the  material  is  categorized  as  orthotropic.  Defining  the  material  properties 
in  these  orthogonal  or  principle  directions,  E1U2=E2212=Q.  This  eliminates 
the  coupling  between  the  normal  stresses  and  the  shear  strain.  Rewriting 
matrix  equation  3.3  with  orthotropic  material  constants,  yields: 


0 11 

®mi 

E\122 

0 

V 

<r22 

-  = 

E2222 

0 

£22 

<V 

sym 

2-^1212. 

en. 

(3.5) 


Equation  3.5  can  be  expressed  in  terms  of  the  vessel  wall’s 
engineering  properties  as: 


0n 

(  1 

v2i&n  0 

V 

022 

.  = 

- 1 - 

E22  0 

^22 

.012  . 

l,  V12V21  J 

sym  2G12(1- v12v21) 

,ei2. 

where  Eu  and  E22  are  the  Young’s  Modulii  and  vI2  and  v2I  are  the  Poisson’s 
Ratios  in  the  1  and  2  coordinate  directions.  These  four  engineering 
constants  are  related  by  the  expression  vI2E22=v2lEn. 

Inverting  and  expanding  matrix  equation  3.6  to  solve  for  the  strains: 


£n 

Cun  Cll22  o 

0n 

£22 

•  = 

Coon  0 

022 

,ei2 . 

sym  2  C1212 

012. 

where  the  compliances  are  defined  as: 


(3.7) 
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(3.8) 


C  -  _ 
'-'1122 


E 


u 


R 


22 


c  - 
^1212  — 


1 

4Gl2 


For  an  isotropic  material,  E11=E33=E,  via=vai=v,  and  Gia=E/[2(l+v)]. 
Thus,  only  two  constants  are  required  to  characterize  the  constitutive 
relationships  for  a  material  such  as  aluminum.  Combining  3.8  and  3.7, 
the  isotropic  longitudinal  strain  can  be  written  as: 


£u  —  g  (On  VG&) 


(3.9) 


and  the  hoop  strain  as: 


C22  “  (°22  V<yil) 


(3.10) 


An  actuator  constructed  from  composites  would,  in  general,  have 
anisotropic  material  characteristics.  Composite  laminates  are 
manufactured  from  individual  plies  which  posses  material  properties 
aligned  in  their  fiber  and  transverse  fiber  directions.  In  order  to  determine 
the  material  properties  in  the  laminate  reference  frame  (1-2),  the 
contribution  from  each  ply  must  be  considered. 

From  Figure  3.4,  a  ply  coordinate  system  (l*-2*)  is  defined  at  an 
angle,  9,  from  the  laminate  coordinate  frame.  In  the  l,-2’  directions,  the 
ply  orthotropic  material  constants,  E' a^fV,  are  defined  as: 


E'nn  = 


Ill. 


E\\22  ~ 


1  _  y'  o' 
X  V,2  V21 

1  —  \t'  \j' 
x  V12  V2l 


E&n  ~~ 


1  ^2^1 


■®1212 


(3.11) 
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FIGURE  3.4:  Typical  Composite  Laminate  Layup 

To  transform  the  ith  ply’s  E'^  into  the  1-2  laminate  frame,  four 
invariant  constants  are  needed  [17,25]. 

/,  =  (1/ 4)[E'mi  +  +2E,/,22‘] 

h  =  a/8)[E,nni+E^02i-2E'n22i  +4  (3 12) 

/3=(l/2  )[E[n;-E^ 
h  =  (1/4)[E1'1I1<  +  E^1  -2E[J  -4 E[J] 

With  these  invariants,  the  ith  ply’s  material  constants  are 
transformed  into  the  1-2  coordinate  system: 

Eutj  =/,  +  I2  +  I3  cos(2  9t)  + I4  cos(4  9;) 

E2Z22  ■/,+/,-/,  cos(2  +  h  cos(4  9t) 

E u22  =h-h  -h  C08(4  V  (3.13) 

Ena  ^h-h  cos(4  9t) 

Eluj  =  (73/2)  sin(2  <Pt)  + 14  sin(4  %) 

E22i2  =  C/s/2)  sin(2  - 14  sin(4  <Pt) 
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With  the  rotated  in  the  laminate  coordinate  frame,  a  thickness 
weighted  stiffness  matrix,  A,  can  be  determined  by  adding  individually  the 
stiffnesses  of  equation  3.13  multiplied  by  the  ith  ply’s  thickness: 


4*.  =  =  'T 

<.-1 


i  -  I 


(3.14) 


The  average  laminate  stresses,  are  defined  as: 


< 

°22 

•=(7) 

.<4. 

v  t  j 

-^liu  ^1122  2A1112 


■^2222  212 


sym 


2A, 


1212  j 


-n 


'22 


'12 


or,  abbreviated: 


(3.15) 


{oL}=|  [Alfel 


(3.16) 


For  a  balanced  laminate,  where  the  number  of  -V  plies  equals  the 
number  of  plies,  AJJJ2=A22I2=0  (i.e.  the  material  is  quasi-orthotropic). 
Inverting  3.16  for  the  strain  vector: 


{£}  =  ([Ar'(aL} 

Assuming  a  balanced  laminate,  equation  3.17  can  be  expanded  as: 


(3.17) 
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(3.18) 


where  d.  —  Aj2j2(A1jj1A2222  ^1122  ^ 

Equation  3.18  is  analogous  to  equation  3.7.  Defining  the  effective 
laminate  compliance  matrix  CL  as: 

[Cl]  =  MA]-1  (3.19) 
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the  effective  or  “smeared”  laminate  engineering  properties  can  be 
determined  from  equation  3.8  as: 


■pL  —  X 
"11  ~  L 

'-'1111 

V12  = 


12 "  4  ch 


1212 


jpL _ 1 

22~cF 


2222 


i/  _  _  pLftL 
k21  “  •c,22'-'1122 


(3.20) 


where  the  elements  of  CL  are  defined  in  equation  3.7. 

With  3.17  and  3.20,  the  longitudinal  and  hoop  strains  are  defined  as: 


en  ~  jgL  (°ii  ^2°^) 


n 


(3.21) 


C22 


(3.22) 


3.3.2  pressure  -Strain  Relationships 

As  illustrated  in  Figure  3.5,  the  longitudinal  and  hoop  stresses 
acting  in  the  walls  of  a  thin  pressurized  cylinder  are: 


(3.23) 


(3.24) 


These  equations  are  derived  by  equilibrating  the  forces  acting  in  the  1  and  2 
directions  (i.e.  the  pressure  force  versus  the  stresses  in  the  cylinder  wall). 

Substituting  3.23  and  3.24  into  the  isotropic  strain  equations  3.9  and 
3.10  yields  the  longitudinal  strain  as  a  function  of  the  internal  pressure: 

PR 

£n  =  ~^U-2v)  (3.25) 
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FIGURE  3.5:  Pressurized  Circular 


and  the  hoop  strain: 


3&3  Derivation  of  the  Actuation  Authority,  Input  Force, 

And  Stroke  Requirements 

From  the  elementary  definition  of  strain,  the  actuation  displacement 
of  the  P-Strut,  SL,  is  given  by: 

SL  =  Ltffen  (3.31) 

where  Ltfr  is  the  effective  length  of  the  actuator.  The  effective  length  is  the 
length  over  which  the  longitudinal  strain  acts  plus  a  contribution  due  to  the 
end  conditions.  Near  the  boundaries  the  pressure-stress  equations  3.23  and 


3.24  are  not  valid.  In  order  to  satisfy  compatibility  at  the  boundaries, 
additional  moment  and  shear  stress-resultants,  acting  out  of  the  1-2 
coordinate  plane,  must  exist.  These  resultants  can  be  solved  using  linear 
Shell  Bending  Theory  which  predicts  a  dramatic  increase  in  en  near  the 
end  conditions  [16,27].  In  essence,  the  clamped-rigid  endcaps  physically 
restrict  the  radial  swelling  or  the  hoop  strain  of  the  vessel  (see  Figure  3.6). 
As  the  hoop  strain  is  diminished  near  the  boundary,  the  local  longitudinal 
strain  is  increased  due  to  Poisson’s  effect. 

Figure  3.7  depicts  the  strains  as  a  function  of  the  position  calculated 
using  a  complete  Shell  Bending  Theory  solution.  The  peaks  in  en  near  the 
ends  of  the  cylinder  correspond  to  the  end  effects  and  the  reduction  in  hoop 
strain  while  the  dominating  plateau  in  etl  is  the  strain  predicted  by  the 
elementary  pressure-strain  equations  3.25  or  3.27.  Therefore,  over  the 
majority  of  the  cylinder  the  simplified  pressure-strain  equations  3.25  and 
3.26  or  3.27  and  3.28  are  valid;  only  near  the  ends  does  the  solution  require 
the  more  complicated  bending  theory  approach.  The  Shell  Bending  Theory 
solution  is  reported  in  Appendix  A  and  analyzed  extensively  (for  isotropic 
and  anisotropic  cylinders)  by  Graves  in  reference  18. 

The  elongation,  dL,  equals  the  longitudinal  strain  integrated  over  the 
length  of  the  P-Strut.  The  peaks  in  eu  at  the  boundaries  are  accounted  for 
by  adding  an  amount  GJto  the  physical  length  of  the  P-Strut,  L,  i.e.: 

Leff  =  L{l+UJ)  (3.32) 
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FIGURE  3.6:  P-Strut  Radial  Swelling  and  End  Effects 


ISOTROPIC  CIRCULAR  CYLINDER 
End  Efforts 


FIGURE  3.7:  Pressurized  Isotropic  Circular  Cylinder  End  Effects 
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3.3.4  volumetric  Changes  and  Input  Force/Stroke 
Requirements 

Since  etl  and  are  both  positive  for  a  positive  pressure  increase,  the 
vessel  grows  with  pressure.  Figure  3.8  depicts,  for  a  circular  cylindrical 
actuator,  the  longitudinal  elongation  or  the  actuation  displacement,  SL.  In 
addition,  the  radial  expansion  or  dilation  is  labeled  SR. 

Th  volume  change  in  the  longitudinal  direction  is  simply  the 
interior  cross-sectional  area  of  the  actuator  multiplied  by  the  actuation 
authority: 

AVl  =  Ae8L  (3.33) 

The  hoop  dilation,  as  drawn  in  Figure  3.7,  equals  the  P-Strut’s 
circumferential  expansion: 

SR  =  Ree„  (3.34) 

The  associated  radial  volume  increase  is  given  by: 

AVr  =  [(Rc  +  SR  f  -  Re2  ]  nL  -  2  kRc8rL  (3.35) 

From  Figure  3.7,  the  bending  theory  end  effect  has  little  influence  on  the 
hoop  strain.  Therefore,  the  hoop  strain  is  assumed  effective  over  a  length  L 
and  not  magnified  nor  reduced  by  the  end  effect. 

The  final  volume  to  consider  is  the  volume  change  of  the  fluid  due  to 
the  internal  pressure.  The  compression  of  the  fluid  volume,  Vf ,  as  a 
function  of  the  pressure  equals: 

4v'-¥-^r 

where  Br  equals  the  bulk  modulus  of  the  fluid. 

Summing  the  volume  changes  together  yields: 
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(3.37) 


AV total  ~  AV I  +  AVg  +  AV f 

As  illustrated  in  Figure  3.9,  to  maintain  a  fluid  balance  within  the  P-Strut, 
the  total  sum  of  the  volume  changes  in  the  cylinder  must  equal  the  volume 
change  in  the  fluid  reservoir.  The  fluid  in  the  bellows  is  assumed 
incompressible  and  the  cross-sectional  area  of  the  bellows,  Ab,  is  assumed 
constant.  (The  relatively  small  volume  of  fluid  in  the  bellows  experiences 
an  insignificant  amount  of  compression.)  Thus,  the  input  force  stroke  is 
defined  as: 

AV, 


X  —  v  total 

b  A 


(3.38) 


The  input  force,  Fb,  is  related  to  the  stiffness  of  the  flexible  bellows 
and  the  internal  bellows  pressure.  Defining  the  spring  stiffness  of  the 
bellows  as  Kb  and  noting  that  the  displacement  equals  Sb,  the  input  force 
required,  without  the  fluid,  is  K^.  With  the  fluid  present,  the  opposing 
pressure  force  equals  PAb.  Summing  these  together  yields  the  required 
input  force: 

Fb  =  Kb8b  +  PAb  (3.39) 

3.3.5  Static  equations  for  an  isotropic  Circular  Cylinder 

Substituting  equation  3.25  into  3.31  and  3.26  into  3.34  yields  for  an 
isotropic  P-Strut  the  axial  displacement  (i.e.  actuation  authority): 

PR 


and  the  radial  dilation: 


s^L«ma-2v) 


*"=B-f§(2-v) 


(3.40) 


(341) 
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FIGURE  3.8:  Longitudinal  and  Radial  Volume  Changes 
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FIGURE  3.9:  Total  Volume  Change  and  Input  Stroke 
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Inserting  equations  3.33,  3.35,  and  3.36  into  3.37  yields: 

AVtotal  =  8LAe  +  8R2nRcL  +  £-AcL  (3.42) 

Bf 

By  substituting  3.40  and  3.41  into  3.42  and  simplifying  using  3.1,  the  total 
volume  change  as  a  function  of  the  pressure  is  determined  as: 

=  ™<{||[M1-2v)+2L(2-  v)]+-|4  &43) 

From  equation  3.38  and  3.43,  the  input  stroke,  8b,  equals: 

*.  =  PVA^[L^1-2v)*2U2-  v)\+±\  '?•  44) 

where  equals  the  ratio  of  the  interior  cross-sectional  area  of  the  cylinder 
to  the  effective  area  of  the  bellows  as  defined  in  equation  3.2.  Substituting 
equation  3.44  into  equation  3.39  yields  the  input  force  as  a  function  of  the 
pressure: 


Fb  =  P 


2  Et 


[4/r(l-2v)  +  2L(2-v)]  +  i^ 


+4, 


(3.45) 


Solving  equation  3.40  for  the  pressure,  P,  yields: 

P  =  SL  2tE 


(3.46) 


LtfrRS  l-2v) 

Noting  the  definition  of  the  isotropic  strain  ratio,  rj  in  equation  3.29 
and  using  equation  3.46,  equation  3.44  can  be  rewritten  as: 


&=$,*#'  1  + 


2  tEL 


2  L 


{l-2v)LtffRcBf 


(3.47) 
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To  determine  the  input  force  as  a  function  of  the  actuation  authority, 
equations  3.46  and  3.47  are  substituted  into  3.39  to  yield: 


2  L  t  2fffL 
+  j  +  (l-2v)LtfrRcBf 


2tE\ 

+  (l-2v)LtfrR/ 


(3  48) 


Therefore,  the  actuation  authority  can  be  determined  as  a  function  of 
the  input  stroke  and  input  force  by  inverting  3.47  and  3.48  respectively,  i.e.: 


,  2  L  2  tE  L 

1+- — —  +  • 


Leffrt  (1-2  v)LttrRcBr 


(3.49) 


and 


Kb^b 

\  2  L  2  tEL 

11  1 

2  tE\ 

Ltfr  r \  (1-2  v)LttfRcB{ 

+  (1-2  v)LtffRc 

(3.50) 


From  studying  the  preceding  equations,  several  important 
observations  can  be  made.  First,  the  static  performance  equations  3.47 
through  3.48  linearly  relate  the  inputs  force,  Fb,  and  stroke,  to  the 
actuation  displacement,  8t.  Second,  the  area  ratio,  f'*,  amplifies  a  small 
input  force  into  a  larger  pressure  force  in  the  cylinder  (equation  3.45)  at  the 
cost  of  increasing  the  input  stroke,  8,,,  per  elongation,  SL  (equation  3.47  or 
3.49).  This  effect  is  analogous  to  a  mechanical  lever. 

In  equation  3.47,  the  first  term  represents  the  axial  response,  the 
second  the  hoop  swelling,  and  the  third  the  fluid  compression.  If,  as  is 
typical,  LILtf  is  approximately  one  and  the  fluid  is  nearly  incompressible, 
the  necessary  bellows  stroke  scales  as  1+2 /Tj.  Since  the  isotropic  strain 
ratio,  rit  is  less  than  one  for  common  isotropic  materials,  the  hoop  strain 
significantly  diminishes  the  static  performance.  For  aluminum,  ]"» 0.25; 
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therefore,  en  is  only  25%  of  e^.  Even  worse,  the  radial  volume  dilation  is 
approximately  an  order  of  magnitude  greater  than  the  longitudinal  volume 
elongation  (l+2/r,.=9).  This  isotropic  volume  effect  translates  directly  into  a 
much  higher  input  stroke,  8b,  per  unit  4-  The  required  input  force,  Fb,  is 
determined  primarily  by  the  last  term  in  equation  3.48  because  the  bellows 
stiffness,  Kb,  is  generally  small.  The  cylinder  material  properties, 
geometry,  and  effective  bellows  area,  Ab,  affect  Fb.  A  larger  Ab,  and  a  larger 
effective  stiffness  E/( l-2v),  or  smaller  cylinder  dimensions  act  to  increase 
the  required  Fb  per  unit  SL.  A  large  Fb  over  a  long  8b  requires  an  input  force 
generator  which  is  capable  of  a  high  force  over  a  long  stroke — a  limiting 
combination. 

The  detrimental  radial  volume  dilation  and  the  higher  effective 
stiffness  (i.e.  E/(  l-2v))  are  results  of  the  isotropic  material  properties.  If  the 
effective  material's  Young’s  Modulii  were  greater  in  the  hoop  direction 
than  the  longitudinal  direction  and  the  Poisson’s  Ratio  was  reduced,  the 
isotropic  disadvantages  could  be  alleviated.  This  revelation  suggests  the 
use  of  tailored  composite  materials  to  achieve  a  cylinder  stiffer  in  the  hoop 
than  in  the  longitudinal  direction. 

3.3.6  The  Static  equations  for  an  Anisotropic  Cylinder 

The  same  process  which  was  followed  to  derive  the  static  equations 
for  an  isotropic  P-Strut  can  be  generalized  to  the  anisotropic  cylinder. 

The  actuation  authority,  SL,  and  the  radial  dilation,  SR,  are  defined 
from  substituting  equations  3.27  into  3.31  and  3.28  into  3.34. 

<3M) 

and 
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(3.52) 


DD 

<5*  =  /?c-f§H2- vf.) 
R  e  2tE&  21 


The  total  volume  change,  AVMai,  was  given  in  3.42.  Substituting  3.51 
and  3.52  into  3.42  and  using  3.38  yields  the  input  stroke,  Sb,  as  a  function  of 
the  pressure,  P: 


Sb  =  PY. 


2 1 


"11  "9 


■'22 


B, 


(3.53) 


Solving  equation  3.51  for  P,  inserting  the  result  into  3.53,  and 
substituting  rc  from  3.30,  the  input  stroke,  Sf,,  as  a  function  of  the 
elongation,  <5L,  equals: 


Sb  =  dLTb  1+ 


2  L  2tE^L 

Le/f  Pc  ( 1  ~  2  v12  )Lt(rRcBf 


(3.54) 


Substituting  3.54  into  3.39,  and  again  using  3.51  solved  for  P,  yields  the 
anisotropic  input  force  to  actuation  authority  relationship. 


F„  =  5, 


KhYh 


1  + 


2  L 


2  tE^L 


Krc  d  2 vy2)Lt/fRcBf 


2  tE^A 


(1-2  v^)LefrRc 


(3.55) 


Equations  3.54  and  3.55  could  be  inverted  to  solve  for  the  elongation  in 
terms  of  the  input  stroke  and  force.  These  relationships  would  be 
analogous  to  equations  3.49  and  3.50  in  the  isotropic  case. 

The  advantage  of  tailoring  the  laminate  properties  is  obvious  from 
equations  3.30,  3.54,  and  3.55.  For  a  given  axial  stiffness,  increasing  the 
hoop  stiffness,  and/or  adjusting  the  Poisson’s  Ratios,  increases  rc.  In 
addition,  tailoring  the  Poisson’s  Ratios  minimizes  the  required  input  force. 
This  implies,  in  comparison  to  the  isotropic  case,  less  force  and  stroke  will 
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be  wasted  on  radial  dilation  and  consequently,  less  force  and  stroke  will  be 
required  for  longitudinal  actuation. 

Intuitively,  the  hoop  stiffness  can  be  increased  by  manufacturing 
composite  P-Struts  with  the  fibers  running  in  the  hoop  direction.  However, 
blindly  increasing  E ^  versus  EnL  without  considering  other  performance 
issues  can  lead  to  less  optimal  actuator  designs.  Chapter  4  investigates  the 
optimal  composite  properties. 


3.4  ANALYSIS  OF  VISCOUS  FLUID-ORIFICE  DAMPING 

In  the  passive  and  active  modes  of  operation,  the  viscous  fluid  is 
forced  in  and  out  of  an  orifice.  As  the  fluid  flows  through  the  orifice  or 
channel,  it  undergoes  a  shearing  action  which  dissipates  energy.  This  loss 
of  energy  is  the  source  of  damping  in  the  P-Strut. 

The  frictional  shear  losses  can  be  estimated  with  analytical  and 
experimental  channel  flow  solutions.  If  the  flow  is  assumed 
incompressible  (i.e.  constant  density),  steady  (i.e.  the  fluid  properties  are 
constant  in  time),  and  fully  developed  (i.e.  the  velocity  profile  is  uniform 
along  the  length  of  the  channel),  the  energy  dissipated  by  shearing  can  be 
related  to  a  mechanical  dashpot  as  depicted  in  Figure  3.10.  The  dashpot 
coefficient,  C,  associates  the  damping  with  the  physical  properties  of  the 
fluid  and  the  dimensions  of  the  orifice.  The  retarding  force  in  the 
dashpot — C  multiplied  by  the  average  fluid  velocity,  V0 — equals  the 
pressure  force  lost,  AP,  due  to  the  fluid  shearing: 

CV0  =  APAv  (3.56) 
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where  Aa  equals  the  cross-sectional  area  of  the  orifice  (A0=tcD02/ 4).  The 
dashpot  representation  is  used  in  both  the  passive  and  active  modeling  of 
the  following  sections. 


3.4.  l  laminar  Hagen -Poiseuille  channel  Flow 

If  the  fluid  flows  in  smooth  layers  and  no  macroscopic  mixing 
occurs,  the  flow  is  laminar.  Figure  3.10  depicts  the  fluid  in  a  two- 
dimensional  circular  orifice  or  channel.  Due  to  the  forces  acting  on  the 
fluid  including  a  constant  pressure  gradient  and  the  shear  stress,  a 
parabolic  velocity  profile  develops.  Summing  these  forces,  the  pressure 
difference  between  locations  1  and  2  can  be  defined  as  a  function  of  the 
orifice  length,  La,  diameter,  Da,  fluid  absolute  viscosity,  p,  and  average 
channel  velocity,  Va,  as: 

AP  —  Px  —  P2—  32~z-~3-  (3.57) 


where  reference  33  or  53  details  the  derivation. 

From  an  energy  balance  the  differential  pressure  can  be  expressed  as 
a  function  of  the  head  loss,  hlam,  and  the  fluid  density,  p,  as: 

AP  =  /zfamp  (3.58) 


From  equations  3.57  and  3.58: 


Kim  ~  32 


ij£L 

A  pvj>. 


=  32 


KYI 

A  Re 


(3.59) 


The  Reynold’s  Number,  Re=pVa/  p,  is  a  non-dimensional  parameter  which 
characterizes  the  flow.  If  Re<  2300,  the  flow  is  considered  laminar. 

^lam  can  be  defined  by  combining  equations  3.56  and  3.57  as: 
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FIGURE  3.10:  A  Circular  Orifice  with  Viscous  Fluid  Flow 
And  the  Mechanical  Dashpot  Representation 
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Cfam  =  8  npL  (3.60) 

Therefore,  Ctam  is  a  function  of  the  length  and  fluid  absolute  viscosity  with 
units  of  force-seconds  per  length.  Interestingly,  the  orifice  diameter  does 
not  influence  Cto(B. 


3.4.2  Turbulent  Fully-Developed  Channel  flow 


The  fluid  is  considered  turbulent,  such  that  macroscopic  mixing  of 
the  fluid  layers  occurs,  if  the  Reynold's  Number  is  greater  than  2300.  By 
virtue  of  the  nature  of  turbulence,  the  analytical  solutions  available  are  for 
very  specific  mixing  problems.  Blasius  proposed  the  following  equation  for 
the  turbulent  head  loss  in  a  smooth  pipe  [33]: 


0,1582  4 
Re 0  25  D0  ° 


(3.61) 


where  i?e<106.  From  the  definition  of  C: 


Cturb  ~ 


0.1582  L 
Re025  D0P  °  ° 


0.0396 

Re025 


nL0D0pVo 


(3.62) 


The  turbulent  dashpot  coefficient,  C^,  is  a  function  of  the  orifice  diameter 
and  the  average  velocity  within  the  turbulent  channel  flow.  CtuH)' s 
dependence  on  V0  makes  applying  CUrb  difficult,  since  the  mechanical 
dashpot  force  would  be  a  non-linear  function  of  Va.  A  non-linear  dashpot 
coefficient  is  difficult  (and  not  practical)  to  incorporate  into  the  linear 
dynamic  models  developed  for  the  passive  and  active  modes  of  operation  in 
the  following  sections. 


3.4.3  Accelerating  and  Oscillating  channel  flows 

The  laminar  and  turbulent  dashpot  coefficients  (equations  3.60  and 
3.62)  were  derived  assuming  a  fully  developed  flow.  However,  time  and 
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distance  along  the  channel  are  required  to  achieve  a  developed  velocity 
profile.  The  entrance  or  development  length,  L for  laminar  flows  equals 
approximately  D0  Re/ 150  [33].  For  turbulent  flows,  ranges  from  25  to  50 
diameters  [33].  Thus  in  a  short  orifice,  the  flow  never  reaches  a  fully 
developed  state.  The  time  to  achieve  a  developed  velocity  profile,  t 
depends  on  the  pressure  flow  condition.  For  the  preceding  constant- 
pressure-gradient  laminar  and  turbulent  flows,  tdeu  can  be  approximated  by 
dividing  the  estimated  Ldtv  by  the  average  channel  flow  velocity,  Va. 

The  actual  pressure  gradient  in  the  P-Strut  orifice  is  dynamic.  The 
external  forces  (i.e.  the  input  force,  Fb ,  or  the  disturbance  force,  Fd, — see 
Figure  1.2)  generate  dynamic  fluid  flow.  If  the  force  instantaneously 
accelerates  the  fluid  flow,  the  time  to  attain  fully  developed  flow,  fAu,  can  be 
estimated  from  the  non-dimensional  parameter,  tdtu\  proposed  by 
Szymanski  and  published  in  reference  53: 

C  =  (3.63) 


where  u  is  defined  as  the  kinematic  viscosity,  n/  p.  When  the  development 
time,  has  progressed  (measured  from  zero  at  the  instant  the  force  was 
applied)  such  that  0.75,  the  fluid  is  essentially  laminar  [53]. 

If  the  pressure  gradient  (i.e.  the  input  or  disturbance  force)  is 
sinusoidally  applied  at  a  frequency  co,  the  fluid  responds  sinusoidally.  Sexl 
and  Unchida  [53]  derived  a  solution  utilizing  a  non-dimensional  frequency, 
co^',  where: 


ca D2 
4v 


(3.64) 
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If  ct)Au*  is  small  (cu^/cl),  the  fluid  velocity  is  essentially  quasi-static  and  the 
laminar,  Hagen-Poiseuille  equations  are  applicable  [53]. 

For  example,  a  DOW  Coming  200  series  silicon  fluid — a  fluid  used  in 
the  P-Strut — has  an  kinematic  viscosity  rated  at  15.5  in2/s  (10K  centistokes) 
[15].  If  the  orifice  diameter  is  assumed  to  be  70  mil  (1.8  mm): 

Re  =  (4.6xlO_a)V0 

**„*=(  12. 6xl03)**„  (3.65) 

^•=  (7.9x10^)^ 

The  channel  velocity,  Vol  would  need  be  0.5x10“  in/s  to  classify  as 
turbulent  flow.  Only  59  microseconds  are  required  to  develop  laminar  flow 
from  the  initiation  of  an  input  or  disturbance  force.  Moreover,  the  dynamic 
forcing  frequency,  co,  would  need  be  on  the  order  of  IK  Hz  to  negate  the  use 
of  a  quasi-laminar  solution.  Therefore,  the  laminar  flow  equations  and  Clam 
seem  appropriate  to  the  P-Strut’s  passive  and  active  analyses. 

3.4.4  Other  Fluid/Orifice  Considerations 

Certainly,  the  P-Strut’s  hand-drilled  orifice  is  not  smooth.  However, 
due  to  the  high  viscosity  of  the  P-Strut  fluid  and  consequently  the  low  Re,  the 
roughness  effects  are  assumed  to  be  negligible  [33]. 

In  addition,  the  orifice  is  not  tapered  and  thus  the  inlet  and  exit  areas 
do  not  uniformly  transition  into  the  reservoir  or  vessel.  As  a  result,  the 
fluid  flow  expands  or  contracts  dramatically  at  the  ends  of  the  orifice. 

These  events  increase  the  fluid  head  loss.  For  a  square-edged  orifice  (i.e. 
flat  exit  with  no  taper)  an  additional  head  loss,  htnd ,  as  great  as  0.25V/  can 
be  present  [33]. 


63 


The  dashpot  equations  developed  in  the  passive  and  active  models  of 
the  following  sections  do  not  consider  the  head  loss  due  to  the  square-edged 
orifice  exits. 


3.5  PASSIVE  PERFORMANCE 


The  passive  performance  of  the  P-Strut  relates  a  collocated  structural 
disturbance,  Fd,  to  an  elongation  (or  compression),  SL.  The  proposed  fluid 
elastic  actuator,  operating  in  its  passive  mode,  exhibits  the  properties  of  a 
viscous  damper.  These  properties  can  be  ascertained  by  examining  the 
mechanical  behavior  and  developing  a  model  of  the  P-Strut. 


3.5.1  An  Analysis  of  the  Mechanical  Behavior 

Figure  3.11  depicts  a  schematic  of  the  P-Strut  subjected  to  a 
disturbance  force,  Fd.  Fd  is  in  equilibrium  with  the  longitudinal  stress  in 
the  cylinder  wall  and  the  pressure  force  inside  the  cylinder: 

Fd^^Aa.-PA  (3.66) 


The  average  hoop  stress  a22L,  associated  with  the  pressure  increase,  was 
reported  in  equation  3.24.  Solving  3.66  for  ouL  and  substituting  ouL  and 
into  equations  3.27  and  3.28  yields  for  the  longitudinal  and  hoop  strains: 


11  2  tEL 


(1-2  vf2)  + 


n 


and,  noting  v12LE22L=v2]LEnL: 


PR, 


£z2  2tEL 


(2-v£)-3F 


22  El  I  Ax 


(3.67) 


(3.68) 
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FIGURE  3.11:  P-Strut  Schematic  in  Passive  Mode  of  Operation 

The  elongation  of  the  cylinder  is  defined  in  equation  3.31  as  the 
longitudinal  strain  multiplied  by  the  length;  Thus,  substituting  3.67  into 
3.31  defines  the  elongation  of  the  P-Strut  as  a  function  of  the  pressure  and 
disturbance  force: 

a'-  =  e"t  =  l§r(1-2vi i)+^r-  <3-69> 

Atlhij  EjjjAqc 

where  L  is  used  instead  of  the  effective  length,  L^,  to  simplify  the 
development.  Typically,  the  effective  length  correction  factor,  CJ  (equation 
3.32),  is  very  small  (i.e.  Leff~L).  Moreover,  if  rigorous,  separate  effective 
length  corrections  would  have  to  be  included  for  radial  swelling  and  axial 
loading. 

The  total  volume  change  for  the  cylinder  due  to  the  combined 
pressure  and  force  loading  can  be  represented  by  inserting  3.67  into  3.33, 
3.68  into  3.35,  and  substituting  the  results  along  with  3.36  into  3.37: 
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total 


=  4 


PRCL 
2  tEfi 


(l-2v*)  +  2 


PRCL 
2  tE^ 


(3.70) 


The  first  three  terms  are  similar  to  the  static  analysis — the  cylinder  grows 
with  increasing  pressure.  The  last  term  indicates  that  a  disturbance  force 
also  acts  to  expand  the  cylinder,  increasing  the  volume. 

Since  within  the  actuator  the  fluid  volume  is  conserved,  the  total 
volume  change  in  the  cylinder  must  be  balanced  by  an  equal  volume  change 
in  the  fluid  reservoir  (see  Figure  3.8).  From  3.38  the  reservoir  stroke,  Sb,  is 
determined  by  dividing  3.70  by  the  effective  area  of  the  bellows,  Ab: 


PRCL 
2  tE 


(l-2vfa)  + 


a 


PRCL 

fl?22 


(2-v£)  +  -F  + 


PL  ,  (1-2 v^)FdL 


EiiK 


(3.71) 


where  equals  the  area  ratio  of  the  cylinder  to  the  bellows  (equation  3.2). 

Assuming  no  input  force  {Fb  =0),  the  pressure  in  the  bellows,  Pb,  must 
be  in  equilibrium  with  the  bellows  static  stiffness,  Kb,  and  stroke,  <%,.  From 
equation  3.39,  Kb§,=-PbAb.  In  addition,  the  pressure  in  the  cylinder,  P, 
must  be  in  equilibrium  with  the  pressure  in  the  bellows  and  the  pressure 
loss  due  to  the  fluid-orifice  interaction.  As  an  equation, 

P  =  Pb  +  AP  =  ~~r^  +  AP  (3.72) 


With  equation  3.56,  equation  3.72  can  be  rewritten  as: 

p  =  k  (3.73) 

4>  Ao 

where  C  (assuming  laminar  flow)  is  defined  in  equation  3.60  and  V0  is  the 
velocity  of  the  fluid  flowing  from  the  cylinder  through  the  orifice  and  into 
the  bellows.  The  velocity  equals  the  rate  of  bellows  displacement  multiplied 
by  the  area  ratio  of  the  bellows  to  the  orifice,  'Pob  (equation  3.2).  The  area 
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ratio  represents  the  velocity  increase  as  the  flow  leaves  a  larger  reservoir 
area  and  enters  a  small  orifice.  In  Laplace  notation,  V0=-lFobs6b. 
Substituting  into  3.73  and  multiplying  by  the  area  of  the  cylinder,  the 
pressure  force,  PAC,  can  be  expressed  as: 

PAC  =  SbAc(&  +  (3.74) 

Rearranging  3.74  and  simplifying  with  the  area  ratios  defined  in  equation 
3.2,  the  bellows  stroke  as  a  function  of  the  pressure  can  be  defined  as: 


('f'bKb  +  'r0'r0bCs) 


(3.75) 


Substituting  equation  3.75  into  3.71  and  solving  for  the  pressure  as  a 
function  of  the  disturbance  force  yields: 


P  =  -F, 


(l-2vf2)L 


lhA 


_M.(1_2v,)  +  M(  2-vM+A 

2 1#/1  2  a)  tEj±{2  ^)+Bf 


'f/b('rbKb  +  'r0'f'obCs) 


(3.76) 


Inserting  equation  3.76  into  3.69  and  dividing  the  numerator  and 
denominator  by  Ac  defines  the  elongation  of  the  P-Strut  as  a  function  of  the 
disturbance  force. 


SL  =  Fd 


L 

Eh  A 


(l-2vf2)flgL(l-2vf2)L 

_ 2  tEhAr  EhA « _ 

(l-2vh)RcL  |  (2-v£)fleL  t  L  | _ 1 

2  tEh\  tEL22Ac  BfAc  'f'b('PbKb  +  'f'0'FobCs) 

(3.77) 


Equation  3.77  represents  the  dynamic  compliance  of  the  system  ( SL 
versus  Fd).  Simplifying,  the  expression  into  a  function  of  stiffnesses  with 
physical  definitions,  yields: 
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*  _ F  i  _ £ _ 

L~  d  Ks  111'  1 

Kl  +  Kr  +  Kf+  V*Kt  +  V!Cs 


(3.78) 


where  KL  equals  the  longitudinal  stiffness  associated  with  the  pressure- 
induced  axial  elongation: 


K  -  K*En 
L  L(1-2v2l2) 


Kr  corresponds  to  the  hoop  or  radial  dilation  due  to  the  pressure: 


*  2L(2-4) 


Kf  is  the  fluid  volumetric  stiffness  which  comes  directly  from  the 
compressibility  of  the  fluid: 


K  - 

Kf-~r 


Ks  represents  the  P-Strut’s  axial  or  structural  stiffness: 

_  EnAM 


(3.79) 


(3.80) 


(3.81) 


(3.82) 


Having  defined  these  stiffnesses  based  on  physical  effects,  the 
performance  equations  developed  with  just  the  mechanical  properties  could 
be  redeveloped  with  these  stiffnesses.  For  example,  the  input  stroke, 
could  be  redefined  from  equation  3.71  as: 


^Ki+£+i]+£ 


(3.83) 
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3.5.2  Discrete  passive  Stiffness  Model 


Inverting  the  compliance  transfer  function,  equation  3.78,  and 
simplifying  yields  the  passive  stiffness  transfer  function,  K^,,  which  is  a 
function  of  the  complex  frequency  s : 


F„  KfKf,  (cs  +  kh+  Kp) 

pau  8l  (K2L-KsKP)(cs  +  kh)  +  K2LKP 


(3.84) 


where: 


KLKRKf 


KLKR  +  KRKf  +  KfKL 


(3.85) 


c=*P02C  and  kb=Fb2Kb. 


Figure  3.12  depicts  a  three  parameter  stiffness  model  representative 
of  the  fluid  elastic  actuator’s  passive  performance.  A  dynamic  disturbance 
force,  Fd,  is  transmitted  directly  through  the  stiffness  labeled  K,  and  in 
parallel  through  the  series  spring  and  dashpot,  K3  and  C,  respectively. 

The  transfer  function  relating  the  structural  disturbance  force,  Fd,  to 
the  axial  expansion,  8L,  of  the  modeled  system  equals: 


(s)  =  K, 


{l+~*-)C,s  +  Ka 

KI _ 

C,s+Ks 


(3.86) 


The  springs  K,  and  K2  and  the  dashpot  C,  do  not  correspond  one-to- 
one  with  the  physical  stiffnesses  and  the  dasphot  given  in  equation  3.84. 
Nevertheless,  the  modeled  springs  and  the  dashpot  can  be  expressed  as 
functions  of  the  physical  stiffnesses,  KL,  Ks,  Kp,  kb,  and  c  as  follows: 
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(3.87) 


Kfjkb  +  KP ) 

1  s(Kl-KpKs)kb+KlKp 
K2K?K2 

Ks  =  (Kl  -  KPKS  )[kb (Kt  -  KPKS )  +  K\KP  ] 
1 


From  equation  3.87  and  Figure  3.12,  the  static  load  path,  Kv  is 
primarily  determined  by  the  structural  stiffness  of  the  cylinder,  Ks.  If  the 
fluid  is  not  present  (i.e.  Kp=0),  Kt  reduces  to  Ks.  In  practice,  the  presence  of 
the  fluid  contributes  a  small  amount  of  additional  stiffness.  The  second 
modeled  spring,  K2  is  a  complex  combination  of  all  of  the  physical 
stiffnesses  present  in  the  system.  The  dashpot  coefficient,  CIt  represents 
the  orifice  damping,  c,  scaled  by  a  complicated  stiffness  factor. 


FIGRUE  3.12:  Passive  Discrete  Stiffness  Model 
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3.5.3  PASSIVE  FREQUENCY  CHARACTERISTICS 


Equation  3.86  expressed  in  pole-zero  form  equals: 


„(s)  =  Kl 


DCp 


(0 

_£ 

G). 


S  +  Q), 
\S  +  0)PJ 


(3.88) 


The  magnitude  of  at  s=0  (the  DC  gain),  is  labeled  K^.  KDCp  equals  Kr 
From  equations  3.86  and  3.88,  the  pole,  eop,  and  zero,  a>zt  are  determined  as: 


Sl 

C; 


Ct{Kt  +  Kt) 


(3.89) 


The  passive  frequency  characteristics  of  the  P-Strut  in  terms  of  the 
physical  stiffnesses,  KL,  Kp,  Ks,  kb,  and  c  are  calculated  from  equations  3.87, 
3.88,  and  3.89.  These  properties  are  reported  in  Table  3.1. 

The  frequency  performance  can  be  illustrated  with  a  Bode  Diagram. 
Assuming  a  sinusoidal  displacement  e***,  can  be  written  as: 


K^iw)  =  8UM{ (0)  e'  *(e>)  (3.90) 

where  M(co)  equals  the  magnitude  and  0(g>)  equals  the  phase  of  K^.  Midi) 
and  0(o>)  are  defined  in  Table  3.1.  Figure  3.13  depicts  the  Bode  magnitude 
and  phase  of  KFnu  (equation  3.90) .  This  figure  was  calculated  using  the 
representative  physical  properties  found  in  Table  3.2. 
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TABLE  3.1:  Passive  Frequency  Characteristics 


CHARACTERISTIC 

DC  Stiffness, 

Stiffness  at  Infinity, 

Frequency  of  Zero,  eo2 

Frequency  of  Pole,  Q)p 
Bandwidth,  Bw 
Maximum  Loss  Factor,  rj* 

Frequency,  of 

Bode  Magnitude,  M(co) 

Bode  Phase,  #co) 


PARAMETERS 


K  Kfik  +  Kp) 

s  kb(Ki-KpKs)  +  K2LKp 


K2K 


i^s  _ 


Kl~KPKs 


—  ^DCpP 


K  +  Kp 

c 


(Kj  -  KpK,,  )kh  +  lCf  Kp 
c(K2L-KPKs ) 


co  —  CO 

wp 


P-1 


<»z4P  = 


KocpP 


1+ 


(0 

c*,2 


]Zr' 


tan 


(0-1) 


co 

(O. 


a  CO 


where 


o>p  (Kl-KpK^  +  KlKp 
P  <ot  {kb  +  KP){K2L-KPKs) 
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PASSIVE  PERFORMANCE 
Magnitude  (Stiffness)  vs  Frequency 


Frequency  (Hz) 


PASSIVE  PERFORMANCE 


Frequency  (Hz) 


FIGURE  3.13:  Passive  Performance  (K  )  Bode  Diagrams: 

Magnitude  and  Phase  Versus  Frequency 
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The  passive  performance  follows  a  classic  pattern.  K} _ begins  with  a 

value  of  at  s=0  and  exhibits  a  region  of  increasing  stiffness  until 
leveling  off  at  K^.  This  rise  in  magnitude  is  accompanied  by  an  associated 
bell-shape  change  in  phase. 

Physically,  the  fluid  flowing  through  the  orifice  or  dashpot  begins  at 
DC  with  no  motion.  The  dashpot  has  an  effective  zero  stiffness,  c$=0.  At 
low  frequencies,  the  fluid  is  pushed  through  the  orifice  by  the  displacement. 
Due  to  fluid  drag  in  the  orifice,  the  phase  of  the  displacement  lags  the  force. 
As  the  frequency  increases,  the  fluid  begins  to  “lock-up”  the  dashpot.  The 
phase  lag  of  the  displacement  begins  to  diminish  as  less  fluid  flows  through 
the  orifice.  At  high  frequencies,  the  fluid  is  unable  to  displace  through  the 
orifice  and  the  dashpot  approaches  an  effective  infinite  stiffness.  At  each 
extreme,  the  frequency  dependence  of  Kn..  is  removed  since  cs  becomes 
“constant”.  The  rise  in  the  magnitude  of  K^,  corresponds  to  the  “locking- 
up”  of  the  dashpot. 

Two  important  parameters  are  the  peak  loss  factor,  rj*,  which 
represents  a  measure  of  the  maximum  damping  achievable,  and  the  peak 
frequency,  of,  which  equals  the  frequency  at  which  77*  occurs.  The  loss 
factor  is  commonly  defined  for  viscous  damping  as: 


7j(ft>)  = 


Energy  Dissipated  /  Cycle 
2  k  Maximum  Energy  Stored  /  Cycle 


=  tan  <t>((o) 


(3.91) 


77  is  proportional  to  the  structural  damping  coefficient,  £  (see  Chapter  2). 
For  LTI  systems  with  a  sinusoidal  disturbance,  tj=2£  [12].  By 
differentiating  3.91  to  locate  the  maximum  of  77(a)),  if  and  of  can  be  defined 
as  reported  in  Table  3.1. 
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To  optimally  damp  a  particular  structural  mode  at  a  frequency  m of 
should  be  made  equal  to  0)^  Furthermore,  the  larger  the  damping 
bandwidth,  the  more  tolerable  the  damping  is  to  uncertainty  in  o)i  and 
potentially  the  more  structural  modes  which  can  be  damped.  The 
bandwidth,  Bw,  is  defined  as  the  difference  between  the  pole  and  zero 
frequencies  (i.e.  approximately  the  frequency  region  under  the  bell-shaped 
phase  plot). 


3.5.4  PARAMETERS  STUDY  OF  THE  PASSIVE  PERFORMANCE 


Table  3.2  details  the  results  of  a  parameteric  study  of  the  passive 
model.  The  reference  values  are  based  on  realistic  physical  dimensions 
and  properties  of  an  isotropic  P-Strut.  These  reference  constants  are  used 
in  the  active  analysis  of  Section  3.6  and  again  in  the  optimization  study  of 
Chapter  4.  The  percent  variations  in  the  DC  stiffness,  the  stiffness  at 
infinity,  the  peak  loss  factor  and  frequency,  and  the  bandwidth  were 
calculated.  The  reference  parameters  were  varied  independently  by  ±25%. 
The  percent  difference  is  defined  by: 


Percent  Difference  =  100% 


( Value  -  Reference  Value ) 
Reference  Value 


(3.92) 


Since  an  increase  in  rj*  and  bandwidth  is  beneficial  to  the  damping 
performance,  the  parameteric  study  suggests  increasing  the  pressure 
stiffness,  Kp.  This  equates  to  increasing  the  amount  of  force  (or  strain 
energy)  which  is  transmitted  into  the  fluid  and  thus  into  the  dashpot.  KP 
can  be  increased  by  increasing  KR  which  suggests — as  did  the  static 
analysis — increasing  the  effective  hoop  stiffness,  E22L  (equation  3.80). 
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TABLE  3.2:  A  Parameteric  Study  of  the  Passive  Performance  Model 


REFERENCE  VALUES* 


Kg  K  K^  Kl  Kg  Kf 

(Klbffin)  (Klbffin)  (lbtfin)  (K  lbffin)  (Klbflin)  (KlJffin) 


47.1  8.4 

10.0  117.8 

13.9  26.2 

C 

(lbf  s  /  in) 

% 

( K  in2/in2) 

v. 

(K  in2/in2) 

0.044 

31.42 

204.08 

KoCp  K<xp 

(K  lbf/in)  (Klbffin) 

rf 

Of  By, 

(Hz)  (Hz) 

47.8  48.5 

6.75xl02 

1.59  2. 16x10 2 

PARAMETER 

Kpcp 

PERCENT  VARIATION 

K^  if 

of 

Kl 

+25% 

-0.5 

-1.0 

-35.1 

0.4 

-34.8 

Kl 

-25% 

1.2 

2.2 

74.3 

-0.6 

73.0 

Ks 

+25% 

25.5 

25.9 

25.7 

0.2 

25.9 

Ks 

-25% 

-25.3 

-25.5 

-25.4 

-0.2 

-25.5 

Kp 

+25% 

0.2 

0.7 

40.8 

11.8 

57.4 

KP 

-25% 

-0.2 

-0.7 

-36.7 

-11.7 

-44.2 

h 

+25% 

0.2 

0.0 

-11.8 

13.4 

0.0 

K 

-25% 

-0.2 

0.0 

15.4 

-13.4 

0.0 

c 

+25% 

0.0 

0.0 

0.0 

-20.0 

-20.0 

c 

-25% 

0.0 

0.0 

0.0 

33.3 

33.3 

*The  Reference  Values  are  based  on  an  Isotropic,  Circular  Cylindrical  P-Strut  with 
the  following  properties: 

E  =9  Msi  (62  GPa),  5^=200  Ksi,  t  *10  mils  (0.25  mm),  Re=0.5  in  (1.27  cm), 

L=  6in  (15.2  cm),  v=0.3,  u  =0.0014  Ibf  s/in2  (9.6  Kcp),  and  La=  1.25  in  (31.8  cm) 


Increasing  the  Bulk  Modulus  of  the  fluid,  Bn  increases  Kn  and  stiffens  Kp 
(equation  3.81). 

In  addition,  the  parameteric  study  indicates  that  decreasing  the 
pressure-elongation  stiffness,  KL,  would  increase  the  peak  loss  factor. 
Physically,  this  implies  an  increase  in  the  proportion  of  the  fluid  volume 
change  due  to  axial  deformation  as  opposed  to  swelling  or  compressing. 
However,  the  model  study  also  suggests  increasing  the  structural  stiffness 
of  the  cylinder,  Ks.  Ks  and  KL  are  related  by  equations  3.79  and  3.82, 

Ks  =(1-2  v12l)Kl.  Increasing  Ks  while  decreasing  KL  can  be  achieved  by 
decreasing  the  Poisson’s  Ratio  v12L.  Reducing  vI2L  was  beneficial  to  the 
static  performance  as  well. 

Finally,  the  parameteric  study  shows  that  increasing  the  peak 
frequency,  of,  is  inversely  proportional  to  the  dashpot  coefficient  such  that 
increasing  c  (i.e.  a  more  viscous  fluid  or  smaller  orifice  dimensions)  moves 
the  Bode  curves  towards  DC  without  effecting  the  maximum  loss  factor. 

3.5.5  Comparison  with  the  D-Strut 

The  passive  dampers  which  were  mentioned  in  Chapter  2  have 
frequency  characteristics  and  equations  similar  to  those  developed  above. 
Most  notably,  the  Honeywell  D-Strut,  as  a  viscous  damper,  can  be 
represented  by  a  stiffness  model  analogous  to  Figure  3.12.  A  schematic  of 
the  D-Strut  is  depicted  in  Figure  3.14.  Anderson,  et.  al.  published  a 
thorough  explanation  and  analysis  of  the  D-Strut  in  reference  2. 

Although  both  the  P-Strut  and  the  D-Strut  can  be  identically  modeled, 
the  physical  interpretations  of  the  stiffness  parameters  and  the  geometry  of 
the  devices  differ  significantly.  The  D-Strut  contains  an  outer  and  an  inner 
solid  tube.  For  the  P-Strut  the  outer  tube  is  the  cylinder  wall,  while  the 
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inner  tube  is  the  fluid-pressure  stiffness,  Kp.  A  second  and  perhaps  more 
significant  difference  is  the  orientation  of  the  load  path  parallel  to  the 
cylinder  wall.  The  orifice  and  the  flexible  fluid  reservoir  are  arranged 
inside  the  outer  tube  of  the  D-Strut  in  line  with  a  diaphragm  and  the  inner 
tube.  Since  for  the  D-Strut  both  the  metal  diaphragm  and  fluid  deform 
together,  the  overall  stiffness  of  this  load  connection  is  difficult  to  formulate 
[35].  As  a  solution,  Anderson  calculated  all  of  the  modeled  stiffnesses  by 
fitting  experimental  data  [2].  In  contrast,  the  simplified  geometry  of  the  P- 
Strut  permits  the  direct  derivation  of  the  physical  stiffness  parameters  and 
their  coupling  to  the  modeled  parameters  as  accomplished  in  the  Section 
3.5.2. 


OUTER  TUBE  ORIFICE  FLUID  IN  BELLOWS 


8l 
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FIGURE  3.14:  A  Schematic  of  the  D-Strut 


3.5.6  Comparison  with  piezoelectric  Resistive-Shunted  Struts 

A  piezoelectric  resistive-shunted  device  exhibits  passive  damping 
qualities  analogous  to  the  P-Strut  and  the  D-Strut.  The  piezo  generates  a 
voltage  in  response  to  a  structurally  induced  strain.  The  voltage  is  then 
dissipated  through  a  resistive  circuit  [20]. 
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Piezo  dampers  can  be  modeled  with  the  spring-dashpot  system  of 
Figure  3.12.  The  parameters  of  these  models  are  complex  functions  of  the 
mechanical  and  piezoelectric  properties  of  the  component  materials  and 
the  passive  circuitry,  in  contrast  to  the  physical  stiffnesses  of  the  P-Strut. 
The  damping  loss  factor,  r j,  is  dependent  upon  the  ratio  of  the  open-circuit 
to  closed-circuit  stiffnesses  for  the  piezoelectric  material  [20]  whereas  in  the 
viscous  dampers  r]  is  related  to  the  ratio  of  to  K^.  The  peak  frequency 
is  set  by  the  selected  shunting  resistor  [20].  For  viscous  dampers,  the  orifice 
dimensions  and  fluid  properties  fix  the  of. 

Aldrich,  et.  al.  has  recently  completed  work  on  a  piezoelectric  passive 
damper  which  utilizes  a  mechanical  lever  to  enhance  the  damping 
performance  [1].  The  lever  mechanism  could  be  modified  to  work  as  an 
actuator.  The  concept  is  analogous  to  the  P-Strut  fluid  lever. 

3.5.7  Revised  Passive  Model 

In  order  to  seal  the  P-Strut  vessel,  endcaps  were  epoxied  to  the  vessel. 
The  epoxy  layer  demonstrated,  during  testing,  viscous  damping  at  low 
frequencies.  The  discrete  stiffness  model  was  revised  to  include  the  epoxy 
dynamics. 

The  updated  model  adds  three  parameters,  K^,  K(p>  and  Cep.  The 
meanings  are  analogous  to  the  springs  previously  defined:  K equals  the 
DC  stiffness  of  the  epoxy,  Cep  is  the  modeled  epoxy  mechanical  dashpot 
coefficient,  and  Kep  represents  the  stiffness  in  series  with  the  epoxy  dashpot. 

Figure  3.15  depicts  the  revised  discrete  stiffness  model.  The  three 
parameter  epoxy  model  is  placed  in  series  with  the  static  stiffness  of  the 
cylinder,  Kr  Again  assuming  a  sinusoidal  displacement,  the  revised 
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damping  model  has  performance  characteristics  similar  to  the  original 
model  with  the  addition  of  a  second  pole  and  zero  set. 

Figure  3.16  depicts  the  magnitude  and  phase  responses  of  the  revised 
passive  damping  transfer  function,  K  r  (dashed  lines)  and  the  original 
model,  if^./solid  lines).  The  three  epoxy  parameters,  K^,  Ktp,  and  C,p, 
used  to  calculate  K^'  are  representative  of  the  measured  values  reported  in 
Chapter  6.  The  two  pole-zero  pairs  of  the  revised  model  (one  pair  per 
dashpot)  have  two  frequency  regions  where  stiffening  and  phase  change 
occurs.  However,  the  location  of  the  first  zero  and  pole,  which  is  inversely 
dependent  on  the  epoxy  dashpot  coefficient,  Cep,  is  at  approximately  the 
same  location  as  that  of  the  second  pole-zero  pair  which  is  due  to  the  orifice 
damping,  c.  Because  of  the  proximity  of  the  poles  and  zeroes  and  because 
the  stiffness  change  associated  with  the  epoxy  is  much  greater  than  that 
associated  with  the  orifice,  the  analytical  epoxy  damping  overwhelms  the 
orifice  damping.  In  short,  the  orifice  damping  is  unobserved. 

Nevertheless,  including  the  epoxy  parameters  with  the  original  three 
modeled  stiffnesses  and  dashpot  does  not  negate  the  trends  established  in 
the  parameteric  study  of  Section  3.5.4.  Rather,  a  change  in  Ctp,  and 
Ktp  have  the  same  effect  on  the  first  zero-pole  pair  as  changing  Cu  Kx  and  K2 
in  the  original  model.  The  first  pole-zero  pair  is  controlled  by  the  epoxy 
parameters  and  the  second  by  the  cylinder  and  orifice  parameters. 
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FIGURE  3.15:  Revised  Passive  Discrete  Stiffness  Model 


REVISED  PASSIVE  PERFORMANCE 
Magnitude  and  Phase  vs  Frequency 
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FIGURE  3.16:  Revised  Passive  Performance  Bode  Plots 
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3.6  DERIVATION  OF  ACTIVE  PERFORMANCE  EQUATIONS 


The  P-Strut  is  a  dynamic  fluid  elastic  actuator.  The  active  mode  of 
operation  relates  a  dynamic  input  force  Fb  or  input  stroke  4  a  resulting 
actuation  displacement  SL  assuming  no  disturbance  force  (Fd= 0). 


3.e.  l  Mechanical  Analysis  of  the  Active  performance 

Figure  3.17  depicts  the  P-Strut  in  its  active  mode  of  operation.  The 
dynamic  input  force,  Fb,  equals  the  pressure  response  in  the  bellows,  Pb, 
and  the  spring  force  of  the  bellows: 

Fb  =  Pb\  +  Kb8b  (3.93) 

The  pressure  inside  the  bellows  and  the  pressure  loss  through  the 
orifice,  AP,  are  in  equilibrium  with  the  pressure  inside  the  cylinder,  P: 

Pb  =  P  +  AP  (3.94) 

where  AP  was  defined  in  equation  3.56. 

The  velocity  in  the  orifice,  Va,  is  designated  as  the  speed  of  the  fluid 
leaving  the  orifice  such  that  a  positive  pressure  in  the  bellows,  Pb,  forces  the 
fluid  into  the  cylinder.  This  velocity  had  the  opposite  sign  of  the  velocity 
used  in  the  passive  modeling  (equation  3.74).  Therefore,  the  Laplace 
definition  of  Va  as  a  function  of  the  bellows  stroke  is  V0=f/o6s4-  Substituting 
3.56  (using  this  definition  of  V0 )  and  3.93  into  3.94,  multiplying  by  the  area  of 
the  cylinder,  Ac ,  and  solving  for  P  yields: 


'f'hFh-8A'f'n'r^Cs+'rhKh 


(3.95) 


The  pressure,  P,  causes  an  elongation,  SL,  which  was  defined  in 
equation  3.51. 


82 


a 


Fb 


FIGURE  3.17:  P-Strut  Schematic  in  Active  Mode  of  Operation 


The  bellows  stroke,  8b,  as  a  function  of  pressure  was  reported 
previously  as  equation  3.53.  Inserting  3.53  into  3.95  and  solving  for  the 
pressure  as  a  function  of  the  input  force,  Fb,  yields: 


P  = 


2 1 


§f-U-2vf,)  +  .?£(2  -vf.) 

"9 


J22 


(3.96) 


Br 


The  elongation.  8L,  as  a  function  of  the  input  force,  Fb,  is  determined 
by  substituting  3.96  into  equation  3.51  and  simplifying.  The  resulting 
equation  gives  the  compliance  or  admittance  of  the  actuator. 


8,  =Fh 


' 

i i/r  2  vf2) 

tfT  2tE^Ac 

' 

l+[¥'?Cs+¥'X]j§ 

L'fr  / 1  .  ov L  )  (■  2L  (2  I'M 
2Vl2)iE^A 

(3.97) 


The  admittance  transfer  function  ( 8JFb )  can  be  written  with  the 
physical  stiffnesses  defined  in  equations  3.79  through  3.82  as: 
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2l 

F„ 


=  A^(s)  = 


EL 


l  +  [cs  +  fc(,]j-Jr  +  -J-  +  -^-| 


K‘l  Kr  Kf 


(3.98) 


where:  K[  =  =  *L-^- 

L*fr  L'fr 

c=  X¥*C  and  kb=%2Kb. 


(3.99) 


The  effective  length,  Lefr,  (see  Section  3.3.3)  is  reintroduced  to  the  axial 
elongation  stiffness  since  the  only  axial  loading  is  due  to  the  pressure  force 
inside  the  cylinder.  In  contrast  to  the  passive  analysis,  using  the  effective 
length  for  the  axial  elongation  stiffness  and  the  physical  length,  L,  for  the 
remaining  stiffnesses  does  not  complicate  the  analysis.  The  terms  in 
equation  3.98  are  otherwise  the  same  as  those  found  in  the  passive  analysis, 
equation  3.78. 

The  active  model  evaluated  at  DC  must  match  the  static  equations 
derived  in  Section  3.3.  Setting  s=0  in  equation  3.98  (which  is  equivalent  to 
ignoring  the  dashpot  loss  factor)  and  rearranging  yields: 


ACTIVE 

DC 


act 


1 

^T*TtKbK’L 

1  1  1  " 

K'L+KR  +  Kf 

(3.100) 


Note  the  strain  ratio  Tc  (equation  3.30),  can  be  written  as  the  ratio  of  the 
radial  dilation  stiffness  (equation  3.79)  to  the  axial  or  elongation  stiffness 
(equation  3.80): 

re  =  iu.  =  -El L  (3.101) 

2  Kl 
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The  static  admittance  was  reported  as  the  inverse  of  equation  3.55. 
Substituting  the  definitions  of  the  physical  stiffness  (equations  3.79  through 
3.82  and  3.99),  equation  3.101  into  3.55,  and  inverting  yields  for  the  DC  static 
admittance: 


STATIC 

Rearranging  the  denominator  in  equation  3.102  yields  the  same  equation  as 
the  DC  admittance  of  the  active  analysis. 

3.6.2  Active  Frequency  Characteristics 

As  illustrated  in  Figure  3.17,  the  input  force,  Fb,  is  non-collocated 
with  the  actuation  displacement,  <5,  .  Due  to  the  absence  of  the  competing 
stiffnesses  which  complicated  the  modeling  of  the  passive  stiffness  transfer 
function  (equation  3.84),  the  active  admittance  transfer  function  \  versus  Fb 
can  be  pictured  with  a  simple  five  parameter  discrete  stiffness  model. 
Figure  3.18  depicts  the  active  spring-dashpot  model. 

The  model  parameters  are  the  physical  stiffnesses  and  dashpot, 
namely  KL,  KR,  Kf  ,  Kb  and  C.  To  rederive  the  admittance,  A^,  from  the 
model,  the  input  force,  Fb  is  related  to  the  input  stroke,  <$, ,  including  the 
proper  area  ratios  for  the  dashpot,  (i.e.  *F0  *f'oA.  from  equation  3.95).  4  is 
written  as  a  function  of  SL  using  the  series  stiffnesses  and  the  ratio  of  the 
bellows  area  to  the  vessel  area,  . 

The  model  depicts  three  parallel  load  paths.  The  first  is  the  bellows 
static  stiffness,  Kb.  The  mechanical  dashpot  represents  the  second  and  the 


frequency  dependent  or  damping  load  path.  The  remaining  load  path 
represents  the  fluid-vessel-pressure  stiffness  and  is  formed  from  KL,  KR, 
and  Kf. 

Reducing  the  model  or  simplify  equation  3.98  yields  the  following 
frequency  dependent  admittance: 


_ w _ 

{kb  +  cs)(KLKR  +  KLKf  +  KRKf )  +  KLKRKf 


(3.103) 


represents  the  amount  of  actuation  displacement  (i.e.  elongation) 
achieved  per  input  force.  The  units  of  A ^  is  length  per  force  or  the  inverse 
of  the  complex  actuation  stiffness. 

The  root  of  the  denominator  of  equation  3.102  represents  the  pole,  o)c, 
of  the  active  performance. 


(oe  = 


kb(KLKR  +  KuKf  +  KRKf )  +  KLKRKf 
c(KLKR  +  KLKr  +  KRKf) 


(3.104) 


If  the  input  force  is  assumed  sinusoidal,  qjc  is  the  comer  frequency  on 
a  Bode  Diagram.  A  sinusoidal  input  force  causes  a  sinusoidal  response  at 
the  same  frequency  but  with  a  different  magnitude  and  a  shifted  phase. 

The  magnitude  and  phase  of  A ^  versus  frequency  is  depicted  in  the  Bode 
plots  of  Figure  3.19.  The  DC  admittance  and  the  l  oll-off  associated  with  the 
pole  are  apparent. 

Physically  the  actuation  force  commands  a  “static”  displacement  at 
low  frequencies  with  magnitude,  A^^.  As  the  frequency  increases  the 
fluid  traveling  from  the  bellows  through  the  orifice  is  impeded  by  the 
viscous  shearing  action.  The  phase  of  the  actuation  displacement  begins  to 
lag  the  input  force  as  the  dashpot  starts  to  “lock-up”.  At  high  frequencies 
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the  dashpot  is  effectively  rigid  with  an  infinite  stiffness,  cs=  <».  Thus,  as 
the  frequency  increases  the  actuation  displacement  rolls-off  to  zero. 

The  bandwidth  is  defined  by  the  frequency  where  the  admittance 
equals  0.707 ADCact  (i.e.  a  -3  dB  reduction  in  the  DC  magnitude).  For  single 
pole  systems,  the  comer  frequency  equals  the  bandwidth  [37].  As  stated  in 
Chapter  2,  to  actively  damp  a  mode  or  to  alter  the  performance  of  a 
structure  at  a  particular  frequency,  the  P-Strut  must  be  capable  of 
commanding  a  displacement  or  force  at  that  frequency.  Therefore,  the 
bandwidth  limits  the  effective  frequency  range  of  the  fluid  elastic  actuator. 

In  the  time  domain  the  bandwidth  dictates  the  speed  of  the  response. 
A  high  bandwidth  (ooc»0)  indicates  a  quick  response  to  an  input  force.  If  a 
step  input  force  is  initiated  and  held,  the  actuation  displacement  will 
respond  by  elongating  the  amount  predicted  in  the  static  analysis.  In  the 
complex  frequency  or  Laplace  domain,  equation  3.103  subjected  to  a  step 
input  can  be  written  as: 


SlM=  n  Nn  Ft(si 

Dy  +  D2S 


N 


unit  step 


Dy  +  D2S 


\s 


(3.105) 


where  N,D„  and  D2  are  apparent  from  the  numerator  and  denominator  of 
3.103.  The  time  to  rise  to  within  5%  of  the  steady  state  or  static  value,  tr ,  can 
be  determined  with  an  Inverse  Laplace  Transform.  The  Inverse  Laplace 
Transform  of  3.105  equals: 

SL(t)  =  £ - ( 1  -  e'%1 )  =  ADCact  ( 1-  e-* )  (3.106) 
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FIGURE  3.18:  Discrete  Active  Performance  (Admittance)  Model 


ACTIVE  PERFORMANCE 
Admittance  Magnitude  and  Phase 
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FIGURE  3.19:  Active  Performance  Admittance  CA^)  Bode  Diagrams: 
Magnitude  and  Phase  Versus  Frequency 
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At  t=0,  the  force  is  initiated  and  the  response  is  zero.  As  time 
increases,  the  static  elongation  A is  approached.  By  definition  of  tr: 

0.95  =  (l-e'"'*)  or  tr=—  (3.107) 

a>c 

Therefore,  increasing  the  bandwidth  or  increasing  the  comer  frequency 
decreases  the  time  to  rise,  tr. 

Nevertheless,  increasing  the  bandwidth  is  not  necessarily  ideal.  The 
benefit  of  the  P-Strut’s  inherent  roll-off  is  that  the  actuator  cannot  feedback 
noise  or  add  energy  to  higher  (and  possibly  unmodeled)  structural  modes. 
This  prevents  the  controller  from  commanding  inputs  through  the  P-Strut 
which  could  destabilize  the  higher  dynamics  of  the  structure. 

3.&3  PARAMETERS  STUDY  OF  THE  ACTIVE  PERFORMANCE 

The  findings  of  a  parameteric  study  on  the  active  performance  model 
are  reported  in  Table  3.3.  For  the  active  model,  the  four  springs,  KL,  KR,  Kf, 
and  Kb,  and  the  dashpot  coefficient,  c,  were  varied  by  ±25%  and  the 
resulting  impact  on  the  DC  admittance  and  the  comer  frequency,  o)c,  or 
bandwidth,  Bw,  is  listed.  The  percent  difference  was  defined  in  equation 
3.90.  The  reference  values  were  calculated  from  the  same  physical 
dimensions  and  properties  used  in  the  passive  parameteric  study. 

Table  3.3  indicates  that  decreasing  the  longitudinal  stiffness,  Klt  and 
increasing  the  radial  dilation  stiffness,  KR,  and  the  fluid  stiffness,  Kf, 
increases  the  DC  admittance  while  only  minimally  effecting  the  pole  and 
zero.  In  the  passive  parameteric  study  of  Section  3.43,  stiffening  KR  and  Kf 
also  increased  the  passive  performance. 
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TABLE  3.3:  A  Parameteric  Study  of  the  Active  Performance  Model 

REFERENCE  VALUES 


Kl  Kr  Kf,  Kf  C 

(Klbflin)  ( K lbf/in)  (lbffin)  (Klbffin)  (lbfs/in) 

117.8  13.8  10.0  26.2  0.044 

*Fb  A-DCact  G>e  0r  ^art 

(K  in2/in2)  (K  in2/in2)  (in/lbf)  (Hz) 

31.42  204.1  1.23x10^  1.59 


PERCENT  VARIATION 

PARAMETER  coe  or 


KL 

+25% 

-19.4 

0.7 

Kl 

-25% 

31.6 

-1.1 

kr 

+25% 

7.0 

6.4 

kr 

-25% 

-9.8 

-7.7 

Kf 

+25% 

3.6 

3.2 

Kf 

-25% 

-5.5 

-4.5 

K 

+25% 

-11.9 

13.5 

h 

-25% 

15.6 

-13.5 

c 

+25% 

0.0 

-20.0 

c 

-25% 

0.0 

33.3 

The  bellows  stiffness,  Kb,  represents  a  direct  path  from  the  input 
force,  Fb,  to  ground,  therefore  decreasing  kb  intuitively  increases  A^.  In 
addition,  kb  strongly  influences  the  bandwidth. 

The  orifice  dashpot  coefficient,  c,  influenced  the  location  of  the  pole, 
coc,  such  that  increasing  c  moves  the  pole  towards  DC  thereby  reducing  the 
bandwidth  of  the  actuator.  The  same  observation  was  true  for  the  passive 
performance. 
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3.7  CLOSING  COMMENTS  ON  THE  PERFORMANCE  EQUATIONS 


The  model  developed  in  Section  3.6  correlated  well  with  the  measured 
data.  Therefore,  the  additional  epoxy  damping,  which  was  discovered 
during  testing  and  discussed  in  Section  3.5.7,  was  not  modeled  in  the  active 
analysis. 

The  three  separate  models  can  be  combined.  The  resulting  model 
would  include  the  elongation  of  the  actuator,  SL ,  and  the  bellows  stroke,  Sb , 
as  a  dynamic  function  of  the  input  force,  Fb ,  and  the  disturbance  force,  Fd. 
The  analysis  and  experimental  work  presented  in  this  thesis  does  not 
require  the  combination  model.  Nevertheless,  the  author  is  planning  to 
present  the  combined  model  in  a  future  publication  as  part  of  a  control 
experiment  involving  the  P-Strut. 
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Chapter  4: 

Optimizing  the  P-Strut  Design 


4.1  CHAPTER  OUTLINE 

In  Chapter  4  the  fluid  elastic  actuator’s  dimensions  and  physical 
properties  are  optimized  with  respect  to  the  static,  passive,  and  active 
performances.  Three  P-Strut  designs  are  considered:  an  isotropic 
aluminum,  a  hybrid  aluminum-composite,  and  a  balanced  composite.  The 
hardware  selected  and  the  final  actuator  configurations  are  detailed. 

4 J2  OPTIMIZING  THE  PHYSICAL  DIMENSIONS 

The  P-Strut  is  primarily  constructed  from  a  circular  cylinder.  For  a 
cylinder,  the  three  physical  dimensions  to  optimize  are  the  radius,  Rc,  the 
length,  L,  and  the  thickness,  t.  Section  3.3.5  reported  the  isotropic,  static 
equations  for  a  circular  cylindrical  P-Strut.  Considering  an  isotropic 
actuator,  Figures  4.1  and  4.2  depict  the  trends  for  Rc  and  L  versus  the  input 
force  required  to  achieve  a  specified 
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INPUT  FORCE  VERSUS  CYLINDER  RADIUS 
30  c  (Isotropic  P-Strut  Static  Performance) 


0  0.5  1  1.5  2  2.5  3 

Radius  (in) 

FIGURE  4.1:  P-Strut  Vessel  Radius  Versus  Input  Force 


INPUT  FORCE  VERSUS  CYLINDER  LENGTH 
12  r  (Lsotropic  P-Strut  Static  Performance) 


Length  (in) 

FIGURE  4.2:  P-Strut  Vessel  Length  Versus  Input  Force 
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elongation  of  0.33  mil.  These  curves  were  calculated  with  the  reference 
constants  listed  in  Table  4.1  while  varying  only  the  parameter  of  interest 
and  the  cylinder  thickness.  Figure  4.1  indicates  that  a  radius  of 
approximately  0,5  in  (1.27  cm)  is  optimal.  The  input  force  versus  elongation 
cost,  Fb  /  SL,  as  represented  in  Figure  4.2,  levels  off  as  the  length  of  the  P- 
Strut  increases. 

Additional  constraints  on  L  and  t  include  the  required  static  stiffness 
(approximately  Ks  from  equation  3.87)  and  the  physical  limits  of  the  optical 
interferometer  testbed,  on  which  the  P-Strut  actuator  is  to  be  applied  [5,36]. 
The  length  of  the  longerons  on  the  interferometer  are  approximately  7  in 
(17.8  cm).  Existing  component  hardware  (i.e.  passive  dampers  and 
actuators)  have  DC  stiffnesses  ranging  from  22.9K  lbf7in  (4  N/pm)  [3]  to 
approximately  85K  lbffin  (15  N/pm)  [5].  Combining  these  considerations 
with  the  definition  of  Ks  (equation  3.82),  the  proposed  P-Strut  designs  have 
an  L  of  no  more  than  7  in  and  a  t  of  4  to  15  mil  (0.10-0.38  mm). 
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VG  THE  MATERIAL  PROPERTIES  WITH  COMPOSITES 


Tailoring  the  material  properties  has  the  potential  *o  enhance  the 
static  performance  by  increasing  the  hoop  stiffness  and  therefore 
decreasing  the  radial  dilation  as  shown  in  Section  3.3.6.  Similarly,  the 
passive  and  active  performances  may  be  improved  by  decreasing  vl2L. 

Table  4.1  lists  the  results  of  adding  successive  graphite-epoxy 
AS4/3501-6  plies  to  an  aluminum  cylinder.  The  wraps  have  the  graphite 
fibers  running  in  the  radial  direction  (i.e.  V  =  90°).  The  table  lists  the 
calculated  laminate  properties,  the  axial  stiffness  of  the  cylinder,  the 


TABLE  4.1:  A  Comparison  of  the  Number  of  Composite  Wraps 
Versus  the  Static  Performance  of  a  Hybrid 
Aluminum-Composite 

Reference  Constants 


En  ~  20.6  Msi  (142  GPa) 

G12  =  0.87  Msi  (6  GPa) 
tcomp  ~  5-28  mil  (0.134  mm)/ply 
E^-9  Msi  (62  GPa) 

L  =  6.0  in  (15.24  cm) 

K„  =  10  lbflin  (1.75  KN/m) 

Ab  -  0.025  in2  (16.1  mm2) 


E&  =  1.42  Msi  (9.82  GPa) 
-  0.3 

=  8  mil  (0  203  mm) 

Yal  =  0.3 
Rc= 0.60  in  (1.52  cm) 
Br=  588  Ksi  (4.1  GPa) 


RESULTS 


Parameters  Number  of  90°  Wraps 


None 

1 

2 

3 

4 

EtlL 

9.00 

6.25 

4.94 

4.18 

3.69 

Msi 

E1 

^22 

9.00 

13.61 

15.60 

16.70 

17.41 

Msi 

Vl2l 

0.30 

0.14 

0.10 

0.08 

0.06 

Ks 

45.5 

52.8 

58.4 

63.8 

69.0 

K  Ibffin 

rc 

0.24 

0.97 

1.51 

2.00 

2.42 

n/4 

7240 

3297 

2885 

2794 

2800 

lbf/in 

4/4 

476 

172 

134 

120 

113 

in/in 

strain  ratio,  rc,  and  the  actuation  costs,  Fb  /  4  and  4/4-  The  fibers  increase 
E22l  and  reduce  vI2L.  Due  to  the  increasing  thickness  of  the  cylinder  wall,  Ks 
is  increased.  After  just  one  ply,  the  strain  ratio,  rc,  is  considerably 
increased  and  the  performance  costs  are  drastically  decreased.  Subsequent 
wraps  improve  the  performance  at  a  significantly  decreasing  rate  of  return 
due  to  the  additional  thickness.  After  the  third  ply,  the  actuation  cost,  Fb  /4> 
increases  as  a  result  of  the  increasing  thickness.  The  “optimal”  hybrid  P- 
Strut  selected  was  an  aluminum  cylinder  with  a  single  90°  composite  wrap. 
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Table  4.2  considers  the  performance  of  several  all-composite  P-Struts. 
To  avoid  normal  stress-shear  coupling  only  balanced  laminates  were 
examined.  To  maintain  a  reasonable  thickness,  the  number  of  plies  was 
limited  to  four — therefore,  layup  symmetry  was  not  possible  with  most 
combinations.  The  results  indicate  that  having  all  fibers  run  in  the  hoop 
direction  have  the  least  actuation  costs,  Fb  and  4/4-  However,  if  all  the 
fibers  run  in  the  hoop  direction,  the  static  stiffness  is  below  the  desired 
range.  The  [90/+60/90]  layup  offers  a  compromise  with  an  acceptable  vessel 
axial  stiffness,  Ks,  as  well  as  low  Fb/4  and  4/4  cost  ratios. 

An  interesting  phenomenon  is  observed  when  attempting  to 
maximize  the  strain  ratio  (equation  3.30).  This  parameter  includes  both 
stiffnesses  and  a  factor  of  laminate  Poisson’s  Ratios.  As  shown  in  Figure 
4.2,  the  optimal  layup  with  respect  to  Fe  has  fibers  slightly  offset  from  the 
hoop  direction.  This  is  due  to  a  complex  Poisson’s  effect  called  “scissoring”. 


TABLE  4.2:  A  Comparison  of  P-Strut  Composite  Layups 
Versus  the  Static  Performance 

RESULTS 


Parameters  Balanced  Composite  Layups 


[90/0/90] 

[90*] 

[90/+75/90] 

[9CVT60/90] 

[90/+45/90] 

EuL 

7.85 

1.42 

1.47 

1.91 

3.51 

Msi 

r  L 
&22 

14.27 

20.60 

18.88 

14.61 

11.92 

Msi 

V  L 

V12 

0.03 

0.02 

0.05 

0.13 

0.19 

79.1 

19.2 

19.9 

25.8 

47.4 

K  Ibf/in 

rc 

0.88 

8.15 

8.53 

5.53 

1.55 

FJSl 

3630 

1070 

1118 

1487 

2997 

Ibf/in 

4/4 

182 

64.2 

64.6 

75.2 

134 

in/in 
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STATIC  PERFORMANCE  COMPARISON 


FIGURE  4.3:  Poisson’s  “Scissoring”  Effect: 

rc  Versus  Composite  Four  Ply  Layups 


4.4  VISCOUS  FLUID  AND  ORIFICE 

The  passive  and  active  analyses  indicated  the  frequency  dependence 
of  the  P-Strut’s  performance  on  the  fluid  and  the  orifice.  The  ideal  fluid 
would  have  an  infinite  bulk  modulus  (i.e.  incompressible)  with  a  high 
viscosity.  The  DOW  Coming  200  Series  silicon  fluids  have  a  high  Bf  and  a 
viscosity  which  can  be  tailored  to  a  desired  specification  [15].  A  less 
expensive  alternative  is  glycerol  which  has  a  higher  Bf  and  a  viscosity  in  the 
range  of  the  selected  silicon  fluids.  Table  4.3  lists  the  properties  of  the  two 
silicon  fluids  and  the  glycerol  fluid  which  was  utilized  in  the  P-Strut  tests. 
These  parameters  span  the  range  of  practical  viscous  fluid  properties  and 
include  the  range  where  most  viscous  oils  are  present.  The  thermal 
expansion  and  variation  of  the  fluid  properties  with  temperature  were  not 
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modeled,  although  the  silicon  and  glycerol  fluids  exhibit  decreasing 
viscosity  with  temperature  [6,15,31]. 

The  orifices  used  in  the  P-Strut  design  ranged  in  diameter  from  30  to 
84  mil  (0.76-2.13  mm)  and  in  length  from  0.8  to  1.5  in  (1.93-3.81  cm).  An 
additional  fixed  fluid  passageway,  inside  the  bellows-orifice  connector,  has 
an  approximate  diameter  of  0.16  mil  (0.41  cm)  and  length  of  0.3  in  (0.76  cm). 
This  connector  orifice  was  included  in  the  calculation  of  the  orifice  dashpot 
coefficient,  C. 


TABLE  4.3:  Properties  of  the  P-Strut’s  Viscous  Fluids 

Parameters  10K  cs  Silicon*  30Kcs  Silicon*  Glycerol**  Units 

Bulk  Modulus,  Br  2.14x10®  2.14x10®  5.88x10®  psiMVol/Vol 

Absolute  Viscosity,  p  1.41xl0'3  4.22X10-3  1.67xl0'3  lbfs/in@R.T. 

Density,  p  0.0352  0.0352  0.0452  lbm/in3 


*  Reference  15.  ’*  Reference  31. 


4.5  HARDWARE  SELECTION 

The  flexible  reservoir  contains  the  fluid  which  flows  through  the 
orifice  during  actuation.  The  flexible  fluid  reservoir  selected  was  a  metal 
bellows.  The  bellows  offered  the  advantages  over  a  piston-type  mechanism 
of  being  simple,  “leak-proof”,  and  essentially  frictionless.  However,  with 
the  bellows  difficulties  arise  in  calculating  the  effective  cross-sectional 
area,  Ab,  since  during  actuation  the  metal  convolutions  of  the  bellows  open 
and  close.  A  Parker  Metal  Bellows  was  selected  which  had  an  inner 
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diameter  of  0.125  in  (3.2  mm),  an  outer  diameter  of  0.375  in  (9.5  mm),  and 
an  expanded  length  of  approximately  0.46  in  (11.7  mm)  [38]. 

In  order  to  achieve  the  input  force,  Fb,  over  the  stroke  required,  a 
DC  solenoid  was  chosen.  The  solenoid  is  an  “off-the-shelf  force  generator 
which  is  capable  of  commanding  reasonable  forces  over  large 
displacements.  Other  actuators  such  as  piezoactuators  had  limited  stroke 
capabilities  and/or  were  not  feasible  due  to  size  considerations. 
Unfortunately,  resistive  heating  in  the  coiled  wire  limits  the  on-time  as  well 
as  influences  the  size,  mass,  and  force  characteristic  of  the  solenoid  ri9.421. 
The  time  delay  associated  with  the  magnetic  flux  in  the  air-gap  limits  the 
frequency. 

The  solenoid  selected  was  a  G.W.  Lisk  DC  push-armature  conical 
pole  face  solenoid  capable  of  a  force  of  41  bf  over  0.2  in  stroke.  (18  N  over  0.5 
cm)  with  a  maximum  on  time  of  4.5  minutes  [19].  The  time  delay  is 
approximately  5  ms,  indicating  the  dynamic  behavior  of  the  solenoid  is  not 
important  in  the  bandwidth  of  the  P-Strut’s  active  performance. 


46  FINAL  P-STRUT  DESIGN 

Figures  4.4  through  4.6  depict  a  three  view  drawing  of  the  P-Strut. 
Figure  4.7  details  the  solenoid-bellows-orifice  connections.  The  length, 
radius,  and  thickness  of  the  three  P-Strut  designs  are  reported  in  Table  4.4. 
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SCALE:  0.6  in  (figure)  to  1  in  (P-Strut) 
All  Dimensions  in  Inches 


—Endcap 

Additional  Length  Required 
For  Instrumentation 
(see  Chapter  5) 

Solenoid  Bolted 
(12-32  Bolt) 
to  Endcap 


Solenoid 


AXIAL  VIEW  POINT 
Figure  4.5 


TOP  VIEW  POINT 
Figure  4.6 


See  Figure  4. 7  for  Detail 
and  Dimensions  of 
Solenoid-Bellows-Orifice 
Connections 


Bellows 

Bellows-Orifice 
Connector  Cylinder 

Orifice 


Bleed  I  Fill  Port 
10-32  Threaded  Screw 

Epoxied  j 

Endcaps  x  ^  (Threaded  3  /  8-24 

Connector) 


FIGURE  4.4:  P-Strut  Side  View 
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Far  Endcap 


6.15 


t 

1.75 


0.373 


SCALE:  0.6  in  (figure)  to  1  in  (P-Strut) 
All  Dimensions  in  Inches 


Solenoid 


Solenoid-Bellows 

Connector 

■  Bellows 

Bellows-Orifice 

Connector 

Orifice 

Vessel 

.  Threaded  Connector 
Bleed /Fill  Ports 


l=k»  © 


0.35 

0.30 


•  Near  Endcap 
Bleed  I  Fill  Ports 

1-0.25 


1.75 


rufo 


.10 


Vessel 


0.20 

0.15 


FIGURE  4.5:  P-Strut  Axial  View  Solenoid 

Endcap 

Solenoid 
Bolted  to 
Endcap 


0.20 

Solenoid 

Attachment 


FIGURE  4.6:  P-Strut  Top  View 


101 


SCALE:  1  in  (figure)  to  1  in  (P-Strut) 
All  Dimensions  in  Inches 


SOLENOID 
Core 

Attachment 
to  Endcap 
12-32  Bolt 
Lock  Nut 

Casing 
Covering 
Coiled  Wire 

Push  Armature 


SOLENOID-BELLOWS 
CONNECTOR 

Cap  Fits  Over. 
Bellows 

BELLOWS 

CONNECTOR 
FLUID  PASSAGEWAY 

BELLOWS-ORIFICE 
CONNECTOR 


ORIFICE 

10-32Nylon 
Threaded  Screw 
(With  Orifice  Drilled 
Through  Center) 
Connects  To  Endcap 


Orifice  Diameter 
(0.030  -  0.084) 


FIGURE  4.7:  Solenoid — Bellows — Orifice  Connections 
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TABLE  4.4:  The  Dimensions  of  the  Three  P-Strut  Designs 
Circular  Cylinder  Dimensions 


Parameters 

Isotropic 

Aluminum 

Hybrid 
Composite- A1 

All-Composite 

[90AS0/90] 

Units 

t  d 

5.9 

7.9 

•— 

mil 

t comp 

— 

10.2 

37.6 

mil 

L 

6.90 

6.65 

5.50 

in 

Ri  (Interior  Radius) 

0.60 

0.60 

0.60 

in 

Rc  (Mean  Radius) 

0.61 

0.61 

0.62 

in 

4.7  “SUB-OPTIMAL”  ASPECTS  OF  THE  P-STRUT  DESIGN 

Certainly,  the  fluid  elastic  actuator  described  and  pictured  above  is 
not  optimal  with  respect  to  mass.  The  endcaps  and  the  attachments  are  not 
optimized  with  regards  to  size  or  shape.  The  DC  solenoid  selected  has  a 
greater  force/stroke  characteristic  and  consequently,  more  coiled  wire  and 
metal  core  than  required.  In  addition,  moving  the  fluid  into  an  annulus 
configuration  (i.e.  inserting  a  stiff  foam  or  secondary  composite  cylinder), 
instead  of  the  fully  filled  cylinder,  would  reduce  the  fluid  mass  while  only 
nominally  effecting  the  overall  performance.  This  concept  is  depicted  in 
Figure  4.8.  Although  a  performance  versus  mass  tradeoff  for  an  annulus 
configured  P-Strut  could  lead  to  a  more  “optimal”  design,  the  basic 
equations  and  the  optimization  studies  conducted  in  this  chapter  would 
remain  relevant  to  the  new  design. 
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Hoop  Stiffened  Cylinder 


FIGURE  4.8:  An  “Optimal”  Vessel  for  the 
P-Strut  Design 
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Chapter  5: 

Manufacturing  and  Experimental 
Procedures 


5.1  CHAPTER  OUTLINE 

Chapter  5  describes  the  processes  used  to  manufacture  the  isotropic, 
hybrid,  and  composite  P-Strut  pressure  vessels.  In  addition,  the 
experimental  objectives,  instrumentation,  and  procedures  are  discussed. 

5w2  MANUFACTURING  PROCESS 

The  key  component  of  a  P-Strut  is  the  pressure  vessel  containing  the 
viscous  fluid. 

5.2.1  Chemical  Milling  Process  for  Thin  Aluminum  Cylinders 

Available  “off-the-shelF  6061-T4  Aluminum  tubing  had  an  outer 
diameter  of  1.25  in  (3.18  cm)  with  a  20  mil  (0.51  mm)  wall  thickness,  t^,.  In 
order  to  reduce  tal,  a  chemical  milling  process  was  utilized. 

A  simple  alkaline  reaction  capable  of  etching  aluminum  was 
published  by  Cotton  in  reference  8  as: 
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2AI  +  2NaOH  +  2H20  ->  2NaA10a  +  3H2 1 
2NoA102  +  4H20  ->  2NaOH  +  Al203  •  3H20 1 


(5.1) 


The  base,  sodiu^  hydroxide,  dissolves  the  aluminum,  releasing  hydrogen 
gas,  and  producing  an  alumina  precipitate.  For  safety,  the  reaction  was 
performed  inside  of  an  air  ducted  hood.  In  addition,  protective  hand  and 
eye  wear  was  worn  since  sodium  hydroxide  is  caustic. 

Due  to  the  exothermic  nature  of  the  reaction,  the  rate  is  highly 
dependent  on  the  amount  of  milling  material  and  solution.  To  determine 
the  proper  concentration  (i.e.  defined  as  mass  of  sodium  hydroxide  per 
mass  of  fluid)  and  the  approximate  time  required  to  thin  the  cylinders  a 
series  of  tests  involving  3  in  (7.6  cm)  cylinder  lengths  was  conducted.  The 
results  are  graphed  in  Figure  5.1.  The  curves  are  exponentially  fit  to  the 
data  which  are  the  expected  reaction  rate  behaviors. 

The  actual  P-Strut  aluminum  cylinders,  which  are  approximately 
6  in  (15.2  cm)  in  length  (see  Table  4.4),  were  milled  in  the  fluid  baths  at 
concentrations  of  approximately  15%  sodium  hydroxide.  The  milling  time 
was  approximately  30  minutes,  which  was  le^s  than  anticipated  based  on 
the  earlier  tests.  The  cylinders  were  removed  from  the  basic  solution  and 
measured  with  a  caliper  frequently.  To  insure  uniform  thinning  the 
sodium  hydroxide  bath  was  kept  agitated  by  hand. 

The  aluminum  cylinder  for  the  isotropic  P-Strut  was  milled  to  an 
average  thickness  of  5.9  mil  (0.15  mm)  with  a  standard  deviation  of  0.4  mil 
(0.01  mm).  The  aluminum  cylinder  for  the  hybrid  P-Strut  was  reduced  to 
7.9  mil  (0.20  mm)  with  a  deviation  of  0.4  mil.  The  thickness  was  measured 
at  the  ends  of  the  cylinders.  Since  attempting  to  measure  the  radius  would 
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distort  the  thin  tube,  the  inner  and  outer  radii  were  estimated  from  their 
pre-etched  values  and  the  post-milled  measured  thickness. 


CHEMICAL  MILLING  TEST  RESULTS 
24  r  For  3  Inch  Aluminum  Cylinder  Teat  Lengths 


Time  (mins) 

FIGURE  5.1:  Chemical  Milling  Sample  Test  Results 


5.2.2  Composite  cylinders  for  the  hybrid  and  All-Composite 

P-Struts 

The  composite  manufacturing  process  was  completed  in  the 
Technology  Laboratory  for  Advanced  Composites  (TEL AC)  at  MIT. 
Reference  28  details  the  materials,  preparation,  and  cure  cycle  followed. 

For  the  hybrid  vessel,  the  milled  aluminum  cylinder  was  hand 
wrapped  in  one  ply  of  AS4/3501-6  35%  resin  graphite-epoxy  composite 
prepreg  with  fibers  in  the  hoop  (90°)  direction.  For  the  all  composite  vessel 
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the  laminates  [90/^60/90]  were  cut  and  laid-up  around  a  non-porous  Teflon 
coated  aluminum  mandrill.  For  both  P-Struts  each  composite  ply  was 
overlapped  at  its  joints  to  allow  some  thermal  expansion  during  curing, 
the  overlaps  are  potential  locations  of  voids  see  Figure  5.2.  Conforming  the 
plies  to  the  cylindrical  shape  increased  the  difficulty  of  the  layup  process. 

The  cylinder  was  covered  with  several  materials  in  preparation  for 
curing.  Figure  5.3  depicts  these  material  layers.  The  hybrid  cylinders 
were  cured  without  a  bleeder  pad.  The  all-composite  cylinder  included  one 
layer  of  bleeder  paper.  A  shrink  fit  material,  which  was  the  outer  layer, 
was  secured  to  the  mandrill  with  flash  tape.  This  layer  contracted  during 
the  cure,  applying  pressure  to  the  composite  cylinders. 

The  cure  cycle  followed  was  1  hour  at  240°F  (116°C)  and  2  hours  at 
350°  (177°C)  with  an  8  hour  post-cure  at  350°F  (177°C).  No  external  pressure 
was  used.  The  oven  temperature  was  set  via  a  dial  and  the  temperature 
was  measured  by  a  thermocouple  inside  the  oven  and  reported  by  a  digital 
readout.  The  ramp  time  is  approximately  5°F  (2.7°C  )  per  minute.  The 
cooling  time  was  not  specified. 

The  cured  tubes  were  removed  from  the  oven  and  the  mandrill  after 
cooling  to  near  room  temperature.  The  composite  cylinder  thickness  was 
measured  at  the  ends  with  a  caliper.  The  dimension  were  reported  in  Table 
4.4.  Since  a  single  composite  laminate  is  approximately  5.3  mil  (0.134  mm) 
thick,  the  measured  composite  thicknesses,  (listed  in  Table  4.5)  were 
significantly  higher  than  anticipated.  The  additional  thickness  was 
assumed  to  be  epoxy  or  possibly  voids.  This  explanation  is  based  on  the 
issue  of  composite  conformability  and  not  having  used  a  vacuum  during  the 
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FIGURE  5.2:  Overlap  of  Composite  Plies  During  Layup 
Cross-Sectional  View 


Flash 

Tape 


ZZ 

m 


23 

■  —  / 


Flash 

Tape 


Ends  Supported 
in  Own 
During  Cure 


Gauranteed 

Non-Porous 

Teflon 


Hollow  1.2in 
Outer  Diameter 
Aluminum 
Mandrill 


Porous  Gauranteed 

Teflon  Non-Porous 

Teflon 


FIGURE  5.3:  Axial  Cross-Section:  P-Strut  Cure  Preparation 
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cures.  To  avoid  exaggerating  the  laminate  stiffnesses,  the  data  analysis 
assumed  the  composites  were  formed  from  plies  of  graphite  fibers  5.3  mil 
thick  surrounded  in  epoxy  tcomp  thick. 

Viewing  the  cured  hybrid  cylinder  from  one  end  indicated  that  the 
cylinder  had  deformed  slightly  into  an  elliptical  shape.  This  was  caused  by 
the  shrink-fit  forcing  the  aluminum-composite  vessel  to  fit  the  slightly 
smaller  mandrill  on  which  the  cylinder  was  cured.  Small  wrinkles  on  the 
inner  and  outer  surfaces  of  the  all-composite  P-Strut  cylinder  were  noted. 

5.23  Attaching  Endcaps 

The  endcaps  pictured  in  Figures  4.4  through  4.6  were  attached  to  the 
cured  cylinders  with  3M  Scotch  Weld  2216  B/A  Gray  Epoxy  Adhesive.  This 
epoxy  has  high  axial  and  shear  stiffnesses  [43].  The  epoxy  cures  at  room 
temperature  in  8-12  hours  with  a  full  cure  in  7  days  [43].  No  testing  was 
conducted  before  waiting  at  least  24  hours. 

5.2.4  Filling  the  Cylinder  wrra  Viscous  Fluid 

Having  sealed  the  P-Strut  vessel,  the  viscous  fluid  (reference  Table 
4.4)  was  forced  into  the  cylinder  under  pressure.  The  setup  used  in 
pictured  in  Figure  5.4.  A  port  on  one  endcap  was  utilized  to  fill  while  a 
second  bled  off  the  air  in  the  cylinder.  In  order  to  remove  the  air  trapped  in 
the  bellows/orifice,  the  bellows  was  filled  such  that  air  would  rise  out  of  the 
bellows  (see  Figure  5.4).  In  addition,  the  bellows  were  pumped  repeatedly 
during  the  filling/bleeding  procedure  in  attempt  to  force  the  air  inside  the 
bellows  out  into  the  cylinder  where  it  could  be  bled  off. 

After  filling  the  cylinder  as  completely  as  possible  with  this 
procedure,  two  screws  were  inserted  into  the  bleed/fill  holes.  The  two 
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screws,  coated  in  Teflon  tape,  were  screwed  into  the  cylinder,  slightly 
pressurizing  the  fluid.  This  was  done  to  reduce  the  size  of  any  residual  air 
bubbles. 


Viscous 

Fluid 


FIGURE  5.4:  Bleed/Fill  Diagram 


&3  INSTRUMENTATION  AND  EXPERIMENTAL  PROCEDURE 


The  objective  of  the  experiments  was  to  demonstrate  the  feasibility  of 
the  P-Strut  as  a  fluid  elastic  actuator  and  to  verify  the  models  developed  in 
Chapter  3. 


Ill 


5.3.1  Instrumentation 

Figure  5.5  illustrates  a  side  view  of  the  P-Strut  with  the  experimental 
instrumentation  labeled.  For  static  testing  six  types  of  data  were  gathered: 
input  stroke,  input  force,  longitudinal  strain,  hoop  strain,  internal 
pressure,  and  elongation.  These  measurements  provided  enough 
information  to  confirm  the  static  analysis  of  Section  3.3.  During  passive 
and  active  testing  only  the  end  displacement  and  the  input  force  (i.e. 
disturbance  force  for  passive  testing,  and  input  force  into  bellows  for  active 
testing)  were  recorded. 

Table  5.1  lists  the  instrumentation  used  in  the  experiments.  The 
data  was  digitally  acquired  and  processed  via  a  2630  Tektronics  Box 
connected  to  a  Gateway  2000  (IBM  compatible)  computer  system. 

The  elongation  was  measured  with  a  Zygote  Laser  Interferometer 
which  is  part  of  a  component  tester.  The  component  tester  is  depicted  in 
Figure  5.6.  The  interferometer  operates  by  sensing  the  difference  in 
distance  between  a  beam  reflected  from  the  end  of  the  P-Strut  and  a 
reference  beam. 

In  the  passive  experiments,  a  Piezo  Instruments  244.30 
Piezotranslator  with  a  load  cell  was  used  to  produce  and  measure  a 
disturbance  force.  In  the  static  and  active  tests,  the  piezotranslator  (i.e. 
disturbance  source)  was  disconnected,  releasing  the  end  of  the  P-Strut  to 
elongate  freely.  The  Tektronics  Box  outputted  a  signal  to  command  the 
solenoid  in  the  active  tests  or  the  piezotranslator  in  the  passive  tests.  These 
signals  were  amplified  by  a  Crown  D- 150 A- II  Amplifier  for  the  solenoid 
and  a  Kepco  BOP-500M  Amplifier  for  the  piezo  actuator. 
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TABLE  5.1:  Instrumentation  List  and  Characteristics 


Measurement  Instrumentation  Input  Sensitivity  Range  Noise  or 

Variable  (Type/Manufacturer)  Voltage  Resolution 


Input  Stroke 

(4) 

Trans-Tek  DC  LVDT  24  V 

20.26V/in 

±lin 

NA 

Input  Faroe 

(Fb) 

PCB  Piezotronics  208 A02  20V 

Loadcells  with  Conditioner 

50.5mV/lb 

±5001b 

±0.0021b 

Dist.  Force 

<f4) 

PCB  Piezotronics  208A02  20V 

Loadcells  with  Conditioner 

50.5mV/lb 

±5001b 

±0.0021b 

Long.  Strain 

(fjl) 

Measurement  Group  Strain  »  3V 
Gages  350ft\  120ft2,  and  120fl3 

2.5V/me 

±0.5me 

±5-20/i£ 

Hoop  Strain 

(£22) 

Measurement  Group  Strain  »  3V 
Gages  350ft1 , 120ft2,  and  120ft3 

2.5V/me 

time 

±5-10/x£ 

Pressure 

(P) 

Omega  PX236  Piezoresistive  10V 
Pressure  Transducer 

lmV/psi 

lOOpsi 

±0.25psi 

Elongation 

(SL) 

SERC  Component  Tester  N  A 

w/  Zygote  Laser  Interferometer 

±0.08mil/V±0.8mil 

±l%/sec 

drift 

*  Strain  Gage  Rosette  on  Isotropic  P-Strut:  CEA-06-125UR-360 

2  Two  Separate  Gages  (Longitudinal/Hoop)  on  Hybrid  P-Strut  EA-06-031DE-120 

3  Strain  Gage  Rosette  on  Composite  P-Strut  CE-06-126UT-120 

(Note:  Different  strain  gages  were  used  due  to  availability  constraints.  The  strain  gagea  were  set  to  the  desired  sensitivity  by  the  gains 
an  the  Vishay  2 120 A  Instrument  and  re-balanced  frequently.  The  ranges  were  set  by  input  voltage  limits  for  the  Tektronka  Baa. ) 


5.3.2  Static  Experimental  Procedure 

The  static  tests  were  conducted  by  inputting  a  voltage  step  into  the 
solenoid  and  observing  the  steady-state  responses  of  the  instruments.  The 
data  was  acquired  approximately  one  second  past  the  time  to  rise  as  defined 
in  Section  3.6.2.  Timing  started  at  the  instant  the  force  was  applied. 

Each  test  collected  data  for  a  single  load  level.  At  least  30  tests  were 
conducted  at  a  variety  of  load  levels  to  characterize  the  static  response.  The 
input  voltage  levels  varied  from  5  VDC  to  20  VDC  which  correspond 
approximately  to  0.2  lbf  (0.9  N)  to  2.5  lbf  (11.1  N).  Data  was  acquired  in  both 
directions  (i.e.  increasing  force  levels  and  decreasing  force  levels). 
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The  solenoid  loadcell  had  an  inherent  AC  coupling  in  the  sensor 
which  meant  the  steady  state  value  decayed  to  zero.  Therefore,  the  input 
force,  Fb ,  (which  was  measured  with  the  loadcell)  was  read  at  the  first 
reasonable  point  after  the  step  input.  Initially  the  loadcell  would  read  a 
spike  in  the  force  due  to  the  impact  nature  of  a  step  input.  However,  the 
spike  would  return  to  a  plateau  before  beginning  to  decay.  This  value  was 
considered  the  correct  measurement  of  the  input  force.  The  laser 
interferometer  system  had  a  minor  (approximately  1%  per  second)  drift 
which  could  have  been  temperature  dependent.  Since  data  was  collected 
within  a  few  seconds,  the  effect  was  considered  negligible. 

5.3.3  Passive  Experimental  procedure 

The  Tektronics  system  generated  a  dynamic  signal  that  was  random 
in  frequency  within  a  specified  bandwidth.  This  signal  was  amplified 
through  the  Kepco  Amplifier.  The  output  drove  a  piezo  actuator  (see  Figure 
5.5  and  5.6)  which  provided  the  disturbance  force.  The  disturbance  force, 

Fd,  was  measured  via  a  loadcell;  the  resulting  elongation  of  the  P-Strut,  8U 
was  measured  by  the  laser  interferometric  system.  Both  measurements 
were  recorded  by  the  Tektronics  box. 

The  Tektronics  box  automatically  acquired  800  points  of  data  over  the 
bandwidth  during  each  measurement  cycle.  The  length  of  time  required 
for  each  cycle  was  defined  by  the  bandwidth  and  the  number  of  data  points. 
The  digital  acquisition  system  calculated,  for  each  cycle,  the  Bode  plot  of 
the  stiffness  transfer  function,  Fd/  8L  (see  Chapter  2  for  definition  of  Bode  plot 
and  reference).  The  recorded  transfer  function  represents  at  least  ten 
cycles  of  data  averaged  automatically  by  the  digital  acquisition  system. 
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5.3.4  Active  Experimental  procedure 

In  the  active  experiments,  a  random  frequency  voltage  signal  was 
generated  by  the  output  of  the  Tektronics  Box,  amplified  by  the  Crown 
amplifier,  and  applied  to  the  solenoid  creating  a  dynamic  input  force,  Fb. 
The  input  force  was  sensed  by  a  load  cell  attached  to  the  solenoid  (see 
Figure  5.5).  The  resulting  dynamic  elongation  of  the  P-Strut  was  measured 
via  the  laser  interferometer. 

The  Tektronics  box  automatically  acquired  at  least  400  points  of  data 
over  the  bandwidth  during  each  measurement  cycle.  The  length  of  time 
required  for  each  cycle  was  defined  by  the  bandwidth  and  the  number  of 
data  points.  The  digital  acquisition  system  calculated,  for  each  cycle  the 
Bode  plot  of  the  admittance  transfer  function,  8JFb .  The  recorded  transfer 
function  represents  at  least  ten  cycles  of  data  averaged  automatically  by  the 
digital  acquisition  system.  The  input  was  constrained  by  the  hardware  to  a 
two  decade  frequency  range.  The  range  was  selected,  based  on  preliminary 
tests,  to  include  the  roll-off  of  the  active  response  in  the  Bode  Diagram. 
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Chapter  6: 

Experimental  Results  and  Correlation 


6.1  CHAPTER  OUTLINE 


In  this  chapter,  experimental  results  for  the  P-Strut  are  presented 
and  correlated.  Three  P-Struts  were  tested  in  the  static,  passive,  and  active 
modes  of  operation. 

The  static  measurements  appear  linear  although  there  is  significant 
scatter  in  some  of  the  data.  The  static  equations  of  Section  3.3  show  good 
agreement  with  the  data.  The  values  of  two  unknown  parameters  were 
ascertained  from  the  data. 

The  passive  data  acquired  indicated  damping  which  was  not  directly 
traceable  to  the  orifice.  The  damping  source  was  modeled  in  Section  3.5  as 
the  epoxy  used  in  attaching  the  endcaps.  The  epoxy  model  parameters  are 
fit  to  the  measured  data  and  the  revised  model  is  discussed. 

The  active  model  correlates  with  the  active  measurements  verifying 
the  modeling  procedure  and  the  analysis  presented  in  Chapter  3.  The 
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modeled  roll-off  of  the  actuator  authority  agrees  with  the  experimental 
results  from  the  isotropic  and  composite  P-Struts. 


&2  STATIC  PERFORMANCE  EVALUATION 

6.2.  l  Static  Experimental  results 

The  static  results  are  depicted  in  Figure  6.1  through  6.14.  The  first 
six  figures  illustrate  the  six  types  of  data  acquired  for  the  isotropic 
aluminum  P-Strut  which  contained  10K  cs  silicon  fluid.  Figures  6.7 
through  6.12  illustrate  the  performance  measures  of  the  hybrid  P-Strut 
filled  with  three  different  fluids.  The  last  pair  of  plots  (Figure  6.13  and  6.14) 
represent  the  static  performance  of  the  composite  P-Strut  filled  with 
glycerol  fluid.  Each  plot  includes  a  line  corresponding  to  the  static  model 
predictions. 

Figure  6.1  depicts  the  measured  longitudinal  strain  versus  hoop 
strain  in  the  isotropic  cylinder.  The  scatter  is  probably  the  result  of  noise  in 
the  longitudinal  strain  gage  noted  in  Table  5.1.  Figures  6.2  to  6.4  show  the 
hoop  strain,  pressure,  and  input  stroke  graphed  against  the  input  force. 
Figure  6.5  illustrates  the  static  performance  of  the  isotropic  P-Strut  as 
measured  by  actuation  displacement  versus  the  input  force.  The  actuation 
displacement  versus  the  input  stroke  is  depicted  in  Figure  6.6.  A  moderate 
input  force  (2  lbf  or  9N)  over  a  moderate  input  stroke  (0.15  in  or  0.38  cm) 
produces  a  small  but  precisely  controlled  elongation  of  approximately  0.5 
mil  (0.13  mm).  Due  to  limitations  in  the  number  of  data  collection  channels 
available  in  the  acquisitional  hardware  not  every  data  type  was  collected 
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during  every  test.  This  accounts  for  the  differing  maximum  elongation 
values  plotted  in  Figure  6.5  and  6.6. 

Figure  6.7  and  6.8  show  the  static  performance  of  the  hybrid  P- Strut 
filled  with  10K  cs  silicon  fluid.  Due  to  the  more  massive  cylinder  wall,  the 
maximum  elongation  is  slightly  less  than  for  the  isotropic  P-Strut.  Figure 
6.9  and  6.10  depict  the  performance  of  the  same  P-Strut  filled  with  30K  cs 
silicon  fluid.  No  significant  change  in  performance  is  noted.  Figure  6.11 
and  6.12  depict  the  performance  of  the  hybrid  P-Strut  filled  with  glycerol.  A 
modest  increase  in  performance  is  noted. 

Figures  6.13  and  6.14  shows  a  dramatic  increase  in  static 
performance  for  the  composite  P-Strut  over  the  aluminum  and  hybrid 
P-Struts. 

6.2.2  Correlation  of  the  Static  Results 

Table  6.1  lists  the  performance  metrics,  elongation  per  unit  force 
(&t/Fb)  and  elongation  per  unit  input  stroke  {8J  $,)•  The  data  values  were 
computed  by  determining  the  slope  of  the  least  squares  fit  of  the  data  plotted 
in  the  figures.  The  model  values  were  calculated  as  follows. 

The  unknown  effective  cross-sectional  area  of  the  bellows,  Ab,  was 
determined  such  that  the  absolute  error  between  the  static  model 
predictions  and  the  least  squared  linear  fits  to  selected  data  was 
minimized.  The  selected  data  consisted  of  the  six  comparisons  of  the  types 
shown  in  Figure  6.1  through  6.6.  An  effective  Ab  was  determined  for  each 
P-Strut  independently. 

Due  to  the  accidental  out  of  specification  distortion  of  the  bellows 
during  testing  of  the  aluminum  P-Strut,  the  effective  cross-sectional  area 
was  higher  than  either  value  determined  for  the  hybrid  or  composite 
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P-Struts.  A  new  bellows  was  used  in  the  composite  P-Strut  which  had  a 
calculated  effective  area  only  slightly  less  that  than  that  of  the  hybrid 
P-Strut.  Due  to  the  uncertainty  in  the  properties  of  the  chemically  milled 
aluminum  tubes,  the  assumed  Young  s  Modulus  for  the  aluminum  tube 
(and  the  aluminum  in  the  hybrid  tube)  was  adjusted  to  be  8.5  Msi  (59  GPa) 
and  a  P  isson’s  Ratio  was  assumed  to  be  0.32.  All  other  material  properties 
and  dimensions  are  found  in  Tables  4.1,  4.3,  and  4.4.  These  values  were 
used  to  compute  the  model  parameters  (from  equations  3.79  through  3.82) 
reported  in  Table  6.2. 

The  agreement  between  the  model  and  data  is  good.  In  most  cases, 
the  model  under  predicts  the  performance  of  the  aluminum  P-Strut  and 
over  predicts  the  performance  of  the  hybrid  P-Struts.  The  model  correctly 
predicts  the  greatly  improved  performance  of  the  composite  P-Strut. 
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Long  Strain  (micro  in/in) 


Aluminum  P- STRUT:  Static  Performance 


Aluminum  P-STRUT:  Static  Performance 


FIGURE  6.5:  Isotropic  Performance:  Elongation  vs.  Input  Force 


Aluminum  P-STRUT:  Static  Performance 


FIGURE  6.6:  Isotropic  Performance:  Elongation  vs.  Input  Stroke 
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Hybrid  Al-Comp  P-STRUT:  Static  Performance 


Input  Force  (Ibf) 


FIGURE  6.7:  Hybrid  Performance  with  10K  cs  Silicon  Fluid: 
Elongation  vs.  Input  Force 

Hybrid  Al-Comp  P-STRUT:  Static  Performance 


Input  Stroke  (in) 

FIGURE  6.8:  Hybrid  Performance  with  10K  cs  Silicon  Fluid: 
Elongation  vs.  Input  Stroke 
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Hybrid  Al-Comp  P- STRUT:  Static  Performance 


FIGURE  6.9:  Hybrid  Performance  with  30K  cs  Silicon  Fluid: 
Elongation  vs.  Input  Force 


Hybrid  Al-Comp  P-STRUT:  Static  Performance 


FIGURE  6.10:  Hybrid  Performance  with  30K  cs  Silicon  Fluid: 
Elongation  vs.  Input  Stroke 
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Hybrid  Al-Comp  P- STRUT:  Static  Performance 


Input  Force  (lbf) 


FIGURE  6.11:  Hybrid  Performance  with  Glycerol  Fluid: 
Elongation  vs.  Input  Force 


Hybrid  Al-Comp  P-STRUT:  Static  Perform  race 


Input  Stroke  (in) 

FIGURE  6. 12:  Hybrid  Performance  with  Glycerol  Fluid: 
Elongation  vs.  Input  Stroke 
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Composite  P-STRUT:  Static  Performance 


FIGURE  6.13:  Composite  Performance:  Elongation  vs.  Input  Force 


Composite  P-STRUT:  Static  Performance 


FIGURE  6.14:  Composite  Performance:  Elongation  vs.  Input  Stroke 
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TABLE  6.1:  Static  Performance  Results 


P-Strut  Cylinders/Fluids 

Performance  Aluminum  Hybrid  Hybrid  Hybrid  Composite 
Metric  10K  ca  Silicon  lOKcs  Silicon  30K  cs  Silicon  Glycerol  Glycerol 


w 

Data 

0.20 

0.15 

0.17 

0.16 

0.28 

tnil/lbf 

w 

Model 

0.14 

0.18 

0.18 

0.22 

0.29 

mil/lbf 

4/4 

Data 

2.95 

2.18 

2.18 

2.70 

4.99 

mil/in 

4/4 

Model 

2.53 

2.81 

2.81 

3.67 

6.99 

mil/in 

TABLE  6.2:  Constants  Used  in  Correlation  of  Static  Performance  Results 

P-Strut  Cylinders/Fluids 

Aluminum 

Hybrid 

Hybrid 

Hybrid 

Composite 

Parameter 

lOKcs  Silicon  lOKcs  Silicon 

30K  cs  Silicon 

Glycerol 

Glycerol 

A, 

0.35 

0.17 

0.17 

0.17 

0.11  in2 

L'ffL 

1.03 

1.01 

1.01 

1.01 

1.01  %  of  L 

13 

13 

13 

13 

13  lbf/in 

kl 

73.7 

67.4 

67.4 

67.4 

58.3  K  lbf/in 

kr 

8.2 

28.8 

28.8 

28.8 

87.2  K  lbf/in 

kf 

33.5 

34.7 

34.7 

95.3 

Ks 

27.7 

46.5 

46.5 

46.5 

50.4  K  lbfi'in 
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03  PASSIVE  PERFORMANCE  EVALUATION 


6.3.1  Passive  Experimental  Results 

Figure  6.15  through  6.18  shows  the  results  of  the  passive  tests  on  the 
three  P-Struts.  The  Bode  magnitude  and  phase  of  the  disturbance  force  per 
unit  axial  elongation  are  shown,  along  with  model  correlations  to  be 
discussed  below.  Figure  6.15  shows  the  passive  performance  of  the 
aluminum  P-Strut  filled  with  10K  cs  silicon  fluid  and  fitted  with  an  84  mil 
diameter  1.1  in  length  orifice.  The  hybrid  P-Strut  containing  10K  cs  silicon 
fluid  with  an  82  mil  diameter  0.9  in  length  orifice  is  pictured  in  Figure  6.16. 
Figure  6.17  depicts  the  hybrid  P-Strut  without  fluid.  The  composite  cylinder 
with  an  82  mil  diameter  0.9  in  length  orifice  and  filled  with  glycerol  is 
depicted  in  Figure  6.18. 

The  behavior  of  the  classic  zero-pole  system  as  discussed  in  Section 
3.4  is  observed  in  all  cases.  Stiffness  at  zero  and  infinite  frequencies,  KOCp 
and  ,  the  peak  frequency,  of,  and  the  maximum  loss  factor,  77*,  are 
listed  in  Table  6.3.  Alarmingly  damping  is  noted  in  the  hybrid  cylinder 
shown  in  Figure  6.17  which  contains  no  fluid. 

6.3.2  CORRELATION  OF  THE  PASSIVE  RESULTS 

To  correlate  the  passive  performance  with  the  measure  data,  the 
effective  stiffness  of  the  component  tester  had  to  be  added  to  the  model 
discussed  in  Section  3.5.  The  connection  of  the  piezotranslator  to  the 
P-Strut  (see  Figures  5.5.  and  5.6)  introduces  a  stiff  in-series  spring  into  the 
load  path  which  decreases  the  measured  passive  performance.  The  spring 
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is  considered  to  have  a  fixed  stiffness  of  225K  lb£/in  (40  N/pm)  which  was 
determined  from  experimental  calibrations. 

Figures  6.15,  6.16,  and  6.18  picture  the  predicted  passive  damping 
due  to  the  orifice.  In  the  isotropic  and  hybrid  cases  the  predicted  damping 
due  to  the  orifice  is  extremely  small  and  barely  noticeably  running  along 
the  frequency  axis.  In  addition,  the  stiffness  is  virtually  unchanged 
throughout  the  frequency  range  (i.e.  the  flat  dotted  line).  In  the  composite 
case  modest  orifice  damping  is  predicted  with  a  peak  near  10  Hz.  The  data 
does  not  agree  with  any  of  the  predictions. 

A  revised  model  was  d'  'loptu  in  Section  3.5.7  assuming  the  epoxy 
used  to  attach  the  endcaps  to  the  cylinder  provided  additional  viscous 
damping.  The  damping  action  of  the  -  poxy  was  assumed  to  be  sensitive  to 
uncontrolled  parameters  such  as  bond  thickness;  therefore,  no  attempt  was 
made  to  calculate  the  parameters  of  the  epoxy  model  from  material 
properties  and  geometry  of  the  bond.  Instead  the  values  of  the  additional 
model  epoxy  parameters  were  determined  by  fitting  the  data.  These  are 
listed  in  Table  6.4.  For  all  cases  except  the  composite  P-Strut,  the  revised 
model  fit  the  data  very  well.  The  composite  results  are  not  well  matched 
even  by  the  revised  model  due  to  the  lack  of  the  predicted  orifice  damping. 

The  passive  results  indicate  that  the  P-Strut  in  its  current 
configuration  is  not  an  effective  passive  damper.  Unintentional  effects, 
possibly  the  epoxy  by  which  the  endcaps  are  attached,  demonstrate 
damping  which  overwhelms  any  orifice  contribution. 
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Phase  (deg)  Stiffness  (K  lbf/in)  Phase  (deg)  Stiffness  (K  Ibf/in) 


Aluminum  P-STRUT:  Passive  Performance 


Frequency  (Hz) 


FIGURE  6.15:  Isotropic  Passive  Performance 


Frequency  (Hz) 


FIGURE  6.16  Hybrid  Passive  Performance 
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Phase  (deg)  Stiffness  (K  Ibf/in)  Phase  (deg)  Stiffness  (K  Ibf/in) 


TABLE  6.3:  Measured  Pasbive  Performance 


P-Strut  Cylinders  /  Fluids 

Aluminum  Hybrid  Hybrid  Composite 


Parameter 

lOKcs  Silicon 

No  Fluid 

10K  cs  Silicon 

Glycerol 

K DCp 

19.91 

19.95 

22.47 

17.35  K  lbffin 

26.14 

27.6 

27.96 

24.67  K  lbf/in 

of 

0.76 

0.68 

0.74 

0.65 

Hz 

Ti 

0.11 

0.15 

009 

0.15 

TABLE  6.4: 

Epoxy  Parameters  Fit  to  Passive  Data 

P-Strut  Cylinders  /  fluids 

Aluminum 

Hybrid 

Hybrid 

Parameter 

10K  cs  Silicon 

No  Fluid 

10K  cs  Silicon 

Glycerol 

K DCtp 

90 

40 

50 

24 

K  lbf/in 

K 

1000 

42 

35 

10 

K  lbf/in 

ctP 

21 

5.6 

5.1 

2.2 

K  lbf  sec/in 
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6.4  ACTIVE  PERFORMANCE  EVALUATION 


6.4.1  Active  Experimental  results 

The  results  of  testing  the  P-Struts  in  the  active  mode  are  shown  in 
Figure  6.19  through  6.21.  The  Bode  magnitude  and  phase  of  the  admittance, 
actuation  displacement  or  axial  elongation  per  unit  applied  force  ( SL/Fb ), 
are  depicted.  Figure  6.20  plots  the  Bode  data  and  model  for  the  aluminum 
P-Strut  with  10K  cs  silicon  fluid  and  an  orifice  with  a  diameter  of  84  mil  and 
a  length  of  1.1  in.  Figure  6.21  depicts  the  data  for  two  hybrid  P-Struts  both 
with  82  mil  0.9  in  length  orifices.  The  left  most  data  corresponds  to  the 
cylinder  filled  with  30K  cs  silicon  fluid  while  the  second  data  set 
corresponds  to  the  same  P-Strut  containing  lower  viscosity  glycerol  fluid. 
Figure  6.21  depicts  the  effect  of  increasing  the  orifice  diameter  (or 
decreasing  the  orifice  damping)  for  the  composite  P-Strut  filled  with 
glycerol  fluid.  The  orifice  diameters  correspond  with  the  order  of  the  roll¬ 
off  data  curves  such  that  the  highest  diameter  orifice  has  the  greatest 
comer  frequency.  Figure  6.21  depicts  the  composite  P-Strut  data  for  five 
orifices:  30  mil  1.5  in,  40  mil  1.1  in, 55  mil  1.1  in,  70  mil  1.1  in,  and  82  mil  0.9 
in  orifices. 

One-pole  roll-off  behavior  is  observed  in  each  case,  except  the  hybrid 
magnitude  data  which  for  unknown  reasons  appears  anomalous.  The 
observed  DC  admittance,  A and  the  comer  frequency  (a)c)  are  shown  in 
Table  6.5.  The  comer  frequency  of  the  data  was  measured  as  the  frequency 
at  which  the  phase  was  45°. 
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The  DC  admittance  (which  corresponds  to  the  static  performance  of 
the  P-Strut)  is  much  greater  for  the  composite  P-Strut  than  for  the 
aluminum.  The  effect  of  the  orifice  size  on  the  roll-off  frequency  is  clearly 
seen  in  Figure  6.21. 

6.4.2  CORRELATION  OF  THE  ACTIVE  RESULTS 

The  predictions  based  on  the  active  performance  equations  developed 
in  Section  3.6.1  are  shown  on  Figures  6.19  and  6.22  for  the  aluminum  and 
composite  P-Struts  respectively.  Figure  6.22  shows  active  performance  data 
of  the  composite  P-Strut  with  different  orifices  and  the  corresponding 
analytical  predictions.  Figure  6.22  depicts  the  composite  P-Strut  data  for 
three  orifices:  30  mil  1.5  in,  40  mil  1.1  in,  and  82  mil  0.9  in  orifice.  The 
parameters  found  on  Table  6.2,  which  includes  the  effective  bellows  area 
calculated  based  on  the  static  results,  were  used  to  compute  the  model 
performance.  No  fitting  of  the  model  parameters  were  done  to  achieve  the 
correlations.  The  agreement  is  excellent  in  all  cases.  This  correlation 
serves  to  verify  the  validity  of  the  active  P-Strut  model. 
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Phase  (deg)  Admittance  (mils/lbf)  Phase  (deg)  Admittance  (mils/lbf) 


FIGURE  6.19:  Isotropic  Active  Performance 


Frequency  (Hz) 

FIGURE  6.20:  Hybrid  Active  Performance  with  30K  cs 
Silicon  and  Glycerol  Fluids 
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Phase  (deg)  Admittance  (mils/lbf)  Phase  (deg)  Admittance  (mils/lbf) 


FIGURE  6.21:  Composite  Active  Performance  with  Increasing 
Orifice  Sizes  (30, 40, 55.  70,  82  mil) 


6.22:  Composite  Active  Performance 


137 


FIGURE 


TABLE  6.5:  Comer  Frequency  and  Active  Performance  Admittance 


P-Strut 

Cylinder 

Fluid 

Orifice 

Length 

K 

in 

Orifice 

Diameter 

D. 

mil 

Corner 

Frequency 

Hz 

DATA  MODEL 

DC 

Admittance 

Adco*« 

mil/lbf 

DATA  MODEL 

Aluminum 

10K  cs  Silicon 

1.1 

84 

1.85 

2.27 

0.13 

0.14 

Hybrid 

30K  cs  Silicon 

0.8 

63 

0.25 

1.05 

0.15 

0.18 

Hybrid 

Glycerol 

0.8 

63 

13.5 

2.77 

0.09 

0.22 

Composite 

Glycerol 

1.5 

30 

0.10 

0.17 

0.23 

0.29 

Composite 

Glycerol 

1.1 

40 

0.55 

0.76 

0.27 

0.29 

Composite 

Glycerol 

1.1 

55 

2.50 

2.69 

0.24 

0.29 

Composite 

Glycerol 

1.1 

70 

4.81 

7.03 

0.27 

0.29 

Composite 

Glycerol 

0.9 

82 

13.7 

15.9 

0.27 

0.29 
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as  DISCUSSION  OF  RESULTS 


The  P-Strut  concept  explained  in  Chapter  1  has  been  demonstrated. 
P-Struts  have  been  tested  in  static,  passive,  and  active  modes. 

The  orifice  damping  is  not  effective  in  the  passive  mode.  With  the 
current  P-Strut  the  expected  effect  is  small  and  it  is  apparently  masked  by 
other  damping  effects — probably  the  epoxy  by  which  the  endcaps  are 
attached.  It  has  been  speculated  that  air  trapped  in  the  cylinder  or  bellows 
could  also  be  responsible  for  lack  of  the  passive  performance.  More 
fundamental  is  the  non-optimal  geometry  of  the  current  P-Strut  for  the 
passive  damping  function. 

The  static  results  confirm  the  linear  nature  of  the  relationships 
between  input  force  and  stroke,  internal  pressure,  cylinder  hoop  and 
longitudinal  strains,  and  axial  elongation.  The  scatter  observed  in  the 
static  data  is  probably  associated  with  data  collecting  methods.  The  data 
acquisition  hardware  is  designed  for  dynamic  experimentation.  The  AC 
loadcell  and  the  laser  interferometeric  component  tester  discussed  in 
Section  5.3.2  cannot  measure  DC  data  for  more  than  a  short  time  before  the 
signal  decays. 

Supporting  this  conclusion,  the  active  test  results  were  very  good. 
The  active  results  verified  the  P-Strut  model  developed  in  Chapter  3.  This 
model  allows  the  geometry  and  material  properties  of  the  P-Strut  to  be  used 
directly  to  calculate  the  responses.  This  capability  aids  in  understanding 
the  physical  principles  behind  the  modes  of  operations  of  the  P-Strut.  This 
knowledge  can  be  used  improve  the  design  or  customize  the  performance. 


The  importance  and  usefulness  of  tailoring  the  material  properties 
and  the  geometry  have  been  demonstrated.  The  optimal  composite  cylinder 
greatly  improved  the  performance  by  constraining  the  undesirable  radial 
deformations.  The  choice  of  fluids  effects  both  the  magnitude  of  the 
performance  via  the  fluid  bulk  modulus  and  the  comer  frequency  in  the 
active  performance  via  the  fluid  viscosity.  The  orifice  diameter  has  a 
demonstrated  strong  effect  on  the  roll-off  frequency. 

The  P-Strut  presented  here,  can  be  used  as  an  actuator  in  a 
controlled  structures  application.  Its  static  stiffness  and  actuation 
authority  are  similar  to  that  of  other  actuators  such  as  a  piezotranslator. 
Advantages  of  the  P-Strut  include  the  linear  response  to  an  input  command 
and  use  of  “off  the  shelf”  hardware.  Conceptually,  the  P-Strut  could  be 
improved  by  the  addition  of  passive  damping,  however  the  geometry  of  the 
fluid  elastic  actuator  would  need  to  be  optimized  for  this  purpose.  For 
example,  the  Honeywell  D-Strut  as  described  in  Section  3.5.5  is  a  fluid 
elastic  mechanism  optimized  for  passive  damping.  The  usefulness  of  the 
P-Strut  in  its  current  configuration  is  to  provide  axial  actuation  within  a 
specified  bandwidth. 

Although  the  bandwidth  limits  the  modes  which  can  be  targeted  for 
damping,  it  beneficially  prevents  the  amplification  of  noise.  Physically,  the 
P-Strut  cannot  command  displacement  or  force  above  a  set  frequency. 
Therefore,  even  if  an  active  control  compensator  is  unstable  at  high 
frequencies  or  the  feedback  noise  is  severe,  the  actuator  cannot  add  energy 
into  the  higher  modes.  Simply,  a  control  system  using  this  fluid  elastic 
actuator  cannot  destabilize  the  structure  at  high  frequencies. 
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Chapter  7: 

Overview  and  Conclusions 


7.1  OVERVIEW 

A  pressurized  elastic  vessel  undergoes  elastic  deformations.  If  the 
pressure  is  controlled,  the  elongation  is  predictable.  This  fluid  elastic 
actuation  concept  was  used  to  develop  the  P-Strut.  Fluid  pressure, 
controlled  by  a  commanded  input  force,  produced  axial  elongation.  The 
static  operation  was  observed  to  be  linear  and  precise.  Dynamic  input  force 
caused  a  dynamic  elongation  which  was  limited  in  bandwidth  due  to  an 
orifice  through  which  the  viscous  fluid  flowed. 

Models  of  the  static,  passive,  and  active  performance  were  developed 
from  mechanical  principles.  The  models  allowed  the  optimization  of  the 
material  and  geometric  properties.  In  addition,  the  models  can  be  extended 
to  the  design  of  other  similar  fluid  elastic  actuator  concepts. 

Three  P-Struts  were  manufactured  and  tested.  The  experimental 
results  verified  the  static  and  active  models.  In  addition,  the  tailored 
composite  properties  of  the  third  P-Strut  resulted  in  dramatically  better 
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performance.  Other  experimental  results  indicated  that  the  frequency 
performance  of  the  fluid  elastic  actuator  can  be  tailored  either  by  the  fluid 
viscosity  or  more  effectively  by  the  diameter  of  the  orifice.  The  passive 
performance  was  masked  by  other  damping  effects  of  uncertain  origin. 
The  analytical  modeling  and  experimental  results  showed  that  the  P-Strut 
geometry  is  not  optimal  for  this  function.  As  an  actuator,  the  P-Strut's 
results  are  excellent. 


7.2  CONCLUSIONS 

The  P-Strut  is  a  functional  fluid  elastic  actuator  with  properties 
comparable  to  other  actuators  over  a  specified  bandwidth.  In  addition,  the 
P-Strut  is  manufactured  from  “off-the-shelf*  components.  The  physical 
principles  by  which  the  P-Strut  operates  has  been  mathematically  modeled. 
These  models  allow  optimization  of  P-Strut  designs. 

The  P-Strut  can  be  effectively  used  in  the  control  of  a  structure.  The 
bandwidth  inhibits  actuation  above  the  frequency  of  the  roll-off  regardless  of 
the  stability  of  the  compensator.  The  linear  nature  of  the  actuator 
simplifies  control  and  allows  the  damping  of  low  frequency  modes.  The  P- 
Strut  can  also  be  used  to  command  and  hold  static  displacements. 
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7.3  FUTURE  WORK 


A  final  demonstration  of  the  usefulness  of  the  current  P-Strut  for 
active  control  is  in  progress  at  SERC  at  MIT.  A  P-Strut  will  be  inserted  into 
the  optical  interferometry  testbed  [5,36],  and  used  to  damp  low  frequency 
disturbances.  The  P-Strut  concept  is  also  being  considered  for  use  in  the 
Stellar  Interferometry  Testbed  Experiment  (SITE)  at  SERC.  For  this 
application,  the  design  will  be  modified  to  optimize  both  the  passive  and 
active  performances.  In  addition,  industry  has  demonstrated  interest  in 
the  fluid  elastic  actuation  concept. 

The  contribution  of  this  work  to  future  fluid  elastic  actuator  is  mostly 
in  the  modeling  of  Chapter  3  and  the  verification  of  the  models  by  the 
results  of  Chapter  6.  The  demonstration  that  the  performance  of  the 
actuator  can  be  tailored  by  properly  selecting  its  components’  material  and 
geometric  properties  is  relevant  to  any  fluid  elastic  actuator  design. 
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APPENDIX  A:  END  EFFECTS  AND  SHELL  BENDING  THEORY 


The  well-known  longitudinal  and  hoop  stress  equations  3.23  and  3.24 
reported  in  Section  3.3.2  for  a  pressurized  circular  cylinder  assume  a 
plane-stress  state.  This  assumption  is  incompatible  with  the  boundary 
conditions.  For  the  P-Strut,  the  end  conditions  are  clamped-clamped. 

Bending  theory  is  capable  of  determining  the  end  effects  and 
incorporates  the  elementary  stress  equations.  However,  the  Bending 
Theory  solution  requires  several  additional  relationships  including  the 
displacements,  strains,  and  forces  in  order  to  acquire  enough  equations  to 
solve  for  the  strains  at  the  boundaries. 

A.  1  The  Boundary  Effects  For  Isotropic  Circular  Cylinders 

A  pressurized,  isotropic,  thin,  closed,  circular  cylinder  with  rigid 
endcaps  is  depicted  in  Figure  A.l.  The  end  conditions  are  assumed  to  be 
clamped-clamped.  Symmetry  and  isotropy  reduce  the  number  of  unknown 
stress-resultants  [16].  After  some  manipulation  [16],  the  key  equilibrium 
condition  is  given  by: 


^  t  O&  p  —  Q 

dx2  Re 

(A.l) 

The  constitutive  laws  are:: 

(A.2) 

E 

°22  “  ^2  _  y2)  +  V 

(A.3) 

Substituting  A.2  into  A.3  yields: 

r,  PRc 

022  =  E  e22  +  v  2  ' 

(A. 4) 

The  bending  moment  stress  resultant,  Af„,  is  defined  as  a  constant 
multiplied  by  the  curvature  of  the  deformed  shell  [16].  The  curvature  of 
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interest  [16]  is  defined  as  the  second  derivative  of  the  radial  displacement, 
w,  which  is  equal  to  negative  the  radial  dilation,  ( w=-SR ).  Therefore, 


M  Et°  d*w-  pd-hhiM 

11  12(1-  v2)  dx2  ~  dx2 


(A.5) 


where  and  by  geometry  £n  =  -  — 

Substituting  A.4  and  A.5  into  A.l,  the  fourth  order  differential 
equation  in  w  is: 

Et  ,v  . 

-aU-TZwHl-v>p=* 


(A.6) 


The  particular  solution  is: 


PR2A  i-~) 


w  =  — 


or  £22  = 


PR'u->’ 


(A.7) 


This  particular  solution  matches  (as  expected)  the  hoop  strain  (equation 
3.26)  determined  from  the  elementary  stress  equations.  In  general,  the 
homogenous  solution  to  a  fourth  differential  equation  is  a  series  of  sines 
and  hyperbolic  sines  which  could  also  be  represented  by  a  complex 
exponential.  The  homogenous  solution,  by  definition,  decays  to  zero. 
Adding  the  particular  and  homogenous  solutions  together  yields. 


where 


w  =  Aj  sin  Xx  sinh  he  +  A,  cos  he  cosh  he - — 

Et 


2PR.H1  -  — )r  •  , 

_  c  2  1  sin  %  c°sh  x  ~  cos  x  sinh  a 


sin  x  cosh  x  -  cos  %  8inh  a 
sin  2x  +  sinh  2x 


(A.8) 


(A.9) 


2PRH1--) r  .  .  . 

-  2  I  cos  x  sinh  x  +  sin  %  cosh  % 


sin  2 x  +  sinh  2x 


(A.  10) 
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(A.ll) 


X  = 


XL  ,  ,<  3U-V2) 

—  and  a  "  R2t2 


From  solving  A.  2  for  the  longitudinal  strain, 


£n  ~ 


PRil-v2) 
2  Et 


-l)  £, 


(A.  12) 


Figure  A.2  plots  these  equations  for  a  cylinder  with  E  =  9  Msi 
(62  GPa),  t  =  10  mils  (0.25  mm),  v  =  0.3,  Re  =  0.605  in  (15.4  mm),  and  L  =  6  in 
(15.2  cm).  The  dramatic  increase  in  the  longitudinal  strain  is  due  to  the 
exponential  terms  or  the  homogenous  part  of  A.  7.  At  the  ends,  the  hoop 
strain  drops  to  zero  as  the  rigid  endcaps  restrain  the  cylinder  from 
swelling.  As  e?2  decreases,  the  Poisson’s  effect  on  the  longitudinal  strain 
(equation  A.  12)  decreases  and  therefore,  en  increases. 

The  total  area  under  the  curve  is  the  actuation  elongation,  SL,  which 
was  defined  in  equation  3.31  as  the  effective  length  multiplied  by  the 
longitudinal  strain.  To  account  for  the  peak  in  eu  near  the  boundaries,  the 
length,  L,  is  increased  by  an  amount  CJ  where  GJis  a  percentage  of  L.  The 
resulting  effective  length,  Ltfn  was  defined  in  equation  in  3.32  and  used  in 
the  analysis  of  Chapter  3. 


CLAMPED 
BOUNDARY 
CONDITION 
(i.e.  Rigid 
Endcaps > 


FIGURE  A.l:  Pressurized  Isotropic  Circular  Cylinder 
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A.2  The  Boundary  Effects  For  Composite  Circular  Cylinders 

If  the  cylinder  is  anisotropic  the  displacements,  the  strains,  and  the 
stress-resultants  are  coupled.  Graves  solved  the  Bending  Theory  problem 
for  a  pressurized,  thin,  clamped-clamped  ended,  balanced-symmetric- 
composite,  circular  cylinder  [7*3].  His  analysis  is  based  on  an  eighth  order 
potential  function  which  efficiently  incorporates  the  displacements, 
strains,  and  stress-resultants.  His  work  included  a  FORTRAN  computer 
program  which  was  reformatted  by  the  author  of  this  thesis  into  MATLAB 
and  verified  against  the  isotropic  analysis  of  Section  A.l. 

The  “optimal”  hybrid,  Graphite/Epoxy-Aluminum,  P-Strut  and  the 
“optimal”  composite — see  Chapters  4  and  5 — are  balanced  circular 
cylinders  but  not  symmetric.  Nevertheless,  if  the  composite  B  matrix — 
which  couples  bending  and  stretching  in  the  material  and  is  zero  for  a 
symmetric  laminate — is  negligible,  Graves  analysis  and  the  program 
remain  accurate  and  consistent  with  the  physics.  The  entries  in  the  B 
matrices  for  the  hybrid  and  composite  P-Struts  are  in  fact  relatively  small. 

For  the  [90AF60/90]  composite  cylinder  proposed  in  Section  4.3,  Figure 
A.3  depicts  the  end  effects.  CJ  for  the  composite  is  less  than  the  isotropic 
cylinder.  Since  the  hoop  stiffness  is  increased  over  the  isotropic  P-Strut, 
is  reduced.  Therefore,  the  peaks  in  en,  which  result  from  the  drop  in  e^2 
through  the  Poisson’s  effect,  are  smaller  in  comparison  to  the  isotropic 
cylinder.  Table  6.2  lists  the  ratio  of  effective  length  to  the  physical  length 
{LJD  for  the  three  P-Struts. 
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FIGURE  A.3:  Hoop  and  Longitudinal  Strain  in  a  Composite  Cylinder 


